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Abstract. The formation of bound states at surfaces of materials with an energy 
gap in the bulk electron spectrum is a well known physical phenomenon. At 
superconductor surfaces, quasiparticlos with energies inside the superconducting 
gap A may be trapped in bound states in quantum wells, formed by total reflection 
against the vacuum and total Andreev reflection against the superconductor. 
Since an electron reflects as a hole and sends a Cooper pair into the 
superconductor, the surface states give rise to resonant transport of quasiparticle 
and Cooper pair currents, and may be observed in tunneling spectra. In 
superconducting junctions, these surface states may hybridize and form bound 
Andreev states, trapped between the superconducting electrodes. In d-wave 
superconductors, the order parameter changes sign under 90° rotation and, as 
a consequence, Andreev reflection may lead to the formation of zero energy 
quasiparticle bound states, midgap states (MGS). The formation of MGS is a 
robust feature of d-wave superconductivity and provides a unified framework 
for many important effects which will be reviewed: large Josephson current, 
low-temperature anomaly of the critical Josephson current, 7r-junction behavior, 
^ TT junction crossover with temperature, zero-bias conductance peaks, 
paramagnetic currents, time reversal symmetry breaking, spontaneous interface 
currents, and resonance features in subgap currents. Taken together these effects, 
when observed in experiments, provide proof for d-wave superconductivity in the 
cuprates. 
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1. Introduction 

Ever since the discovery of high-Tc superconductivity (Bednorz and MiiUer, 1986) 
spin-fluctuation mediated d-wave singlet pairing has been a major candidate for 
superconductive pairing in the cupratcs (Dagotto, 1994; Scalapino, 1995). Ahcady 
from the start it was reahzed that the layered cuprates are low dimensional 
systems with antiferromagnetic undoped parent compounds, with strong electron 
correlations and with highly unusual normal electron properties. Nevertheless, 
phonon-mediated strong-coupling s-wave singlet pairing could not be ruled out because 
of the experimental difficulty to identify essential d-wave signatures which cannot 
also be explained within an s-wavc picture. Looking back, it seems that the 
essential quantitative body of facts, not directly addressing the qualitative question 
of symmetry, was known already by 1990. In particular, one early experimental result 
may be regarded as a clear qualitative signature of d-wave pairing: the so called 
zero-bias conductance peak (ZBCP) observed in the I-V characteristics of electron 
tunneling into thin films of YBaCaCuO already by Geerk et al. (1988). However, the 
d-wave connection was not realized at the time; it was later proposed by Hu (1994) and 
Kashiwaya et al. (1995) and then experimentally verified by Covington et al. (1997) 
in a study of the magnetic field dependence of the splitting of the ZBCP. 

A large effort has gone into finding the most significant evidence for d-wave 
pairing symmetry of the order parameter in the high-Tc cuprates [see e.g. the 
review papers by Dagotto (1994), Scalapino (1995) and Annett et al. (1996) and 
references therein]. The problem is nontrivial because it is necessary to distinguish 
d-wave symmetry from anisotropic s-wave symmetry. A d-wave order parameter 
has nodes at four mutually perpendicular directions of the electron momentum. 
This provides gapless superconductivity in the direction of the nodes, with normal 
excitations existing at all energies, having far-reaching consequences for the materials 
properties of the superconducting state, e.g. linear heat capacity at low temperature, 
large subgap tunnel current, smeared gap structures and high sensitivity to non- 
magnetic impurities. These properties have indeed been found in experiments, 
from the beginning. Moreover, anisotropic gap structures with deep minima in the 
node directions have been established in angle resolved photoemission experiments. 
The problem is, however, that these experimental facts are also consistent with an 
anisotropic s-wave order parameter [see the review by Annett et al. (1996) for an 
analysis of experimc^ntal results] . 

To make an "experimentum crucius" one has to establish that the order parameter 
changes sign when the momentum rotates by 90 degrees. Such experiments have 
earlier been suggested by Geshkenbein et al. (1987) in connection with heavy 
fermion superconductivity, and first realized in tunneling experiments on high-Tc 
superconductors by WoUman et al. (1993). The idea behind these experiments is 
based on the phase sensitive properties of Josephson junctions involving d-wave and 
s-wave, or two misorientcd d-wave, superconducting electrodes. When the junction 
electrodes have different signs of the order parameter, this is equivalent to adding a 
TT-shift to the superconducting phase difference across the junction, and therefore the 
stable state of the junction corresponds to a superconducting phase difference equal 
to ^ = TT (TT-junctions), instead of = like in conventional junctions. 

The conceptually simplest way to create a 7r-junction is to attach a 
superconducting loop to two adjacent faces at right angles of a high-Tc single 
crystal, creating a SQUID circuit with one 7r-junction (Mathai et al., 1995; Wollman 
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et al., 1995; van Harlingen, 1995; Schulz et al., 2000). Using tri-crystal substrates 
one can create a superconducting ring consisting of an odd number of 7r-junctions 
which is able to trap half-integer flux quanta. This was studied by Tsuei et al. 
(1994; 1996; 1997) and Kirtley et al. (19956; 1996) who used a scanning SQUID to 
detect spontaneous half-integer flux quanta in the 7r-junction loop but no flux quanta 
in rings on bi-crystal interfaces; for an overview of the subject, see the recent review 
article by Tsuei and Kirtley (2000a). This is in principle the definitive experiment, 
but it is desirable that it can be repeated by other groups. 

There is however another type of experiments which may provide direct evidence 
for rf-wave symmetry of the order parameter. These experiments exploit the property 
of the d-wave superconductors to form superconducting bound surface states at the 
Fermi energy, so called midgap states, MGS (Hu, 1994). MGS is a direct consequence 
of the angle-dependent sign change of the order parameter, and it cannot exist in s- 
wavc superconductors. Transport measurements which detect the MGS resonance at 
the Fermi energy for particular orientations of the plane of the junction with respect 
to the crystallographic axes will provide crucial evidence for rf-wave superconductivity. 

The problem of MGS is part of a general problem of superconducting surface 
and interface states and their relation to the Andreev states, which determine the 
transport properties of superconducting junctions. This is the topic of the present 
review. 

2. Superconducting surface states 

As is well known, at semiconductor and metal surfaces quasiparticle surface states may 
exist in the band gaps of the bulk material (see e.g. Ashcroft and Mermin, 1976). The 
situation is similar at a superconductor surface: the important difference is, however, 
that on a superconductor surface, electrons are reflected as holes and holes are reflected 
as electrons - Andreev reflection (Andreev, 1964) - while on a normal surface electrons 
are reflected in the normal way as electrons. 

Superconducting quasiparticle surface states were first described by de Gennes 
and Saint-James (1963) for the INS quantum well. More recently it was found that 
the superconducting surface states do not need a normal region at the surface and 
they exist at an arbitrary interface between an s-wave superconductor and an insulator 
(Wendin and Shumeiko, 1996a). The energy of the surface states lies within the energy 
gap of the superconductor and it depends on the properties of the surface. In d-wave 
superconductors, due to the specific symmetry of the order parameter, these levels 
may occur exactly at the Fermi energy (Hu, 1994; Hu, 1998). 

In superconducting junctions, the overlap of the wave functions of the surface 
states builds up stationary states at certain energies, the Andreev bound states. These 
states are thus states of the entire junction, similar to bonding and antibonding states 
in a molecule or in a system of two closely separated metal siirfaces. The name of 
these bound states are associated with Andreev, who first studied them in perfect 
SNS junctions (Andreev, 1965). Due to their ability to carry current, the Andreev 
boimd states play an essential role in dc Josephson transport (Kulik, 1969; Bardeen 
and Johnson, 1972). Moreover, the surface states in superconductor- normal metal 
and superconductor-superconductor junctions introduce resonances in the current 
transport under applied voltage bias. 
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2.1. Andreev reflection 

An important aspect of the superconducting surface states is that they are built up 
from a combination of electron and hole wave functions. This mixing is induced by 
the specific property of the superconductor interface known as Andreev reflection. 

The problem was originally raised by the observation of an anomalously large 
thermal resistance of superconductors in the intermediate state (Mendelssohn and 
Olsen, 1950; Zavaritsky, 1960). It was suggested that the enhancement of the 
thermal resistance was due to some strong additional quasiparticle scattering at the 
NS interfaces (Hulm, 1953; Striissler and Wyder, 1963). However, no microscopic 
mechanism for such a reflection was proposed because the spatial variation of the 
order parameter in the intermediate state is smooth on the scale of the Fermi wave 
length, and it was not clear how such a " soft" potential could provide efficient electron 
reflection. The problem was solved by Andreev who showed that the reflection at NS 
interfaces is due to electron-hole conversion (Andreev, 1964). 

Let us consider a contact between a superconductor and a normal metal, and 
assume for simplicity that the normal electron properties of both electrodes are 
identical. This means that normal electrons will not be reflected by the contact - 
perfect transmission in the normal state. However, electrons in the normal electrode 
which approach the interface with an energy E inside the superconducting gap, 
\E\ < A, cannot penetrate into the superconducting electrode. At the same time, 
there is no normal scattering process at the interface which is able to reverse the 
electron momentum. The problem is solved by converting the incoming electron into 
a reflected hole moving with the same momentum in the opposite direction: since the 
hole energy is given by E'^ = Ep — p^/2m, the group velocity is = dE^ /dp « —vp- 
Such a conversion, Andreev reflection, occurs due to electron-hole correlations within 
the superconducting electrode. 

A straightforward quantum mechanical analysis (see the appendix) shows that 
Andreev reflection is accompanied by a phase shift —j{E) =F x, where the upper 
sign corresponds to the conversion of an electron into a hole, and the lower sign 
to the inverse process. In this expression, jiE) = arccos(i<^/| A|) is the Andreev 
reflection phase shift and x is the phase of the superconducting order parameter. 
The wave functions of the electrons and holes penetrate into the superconductor 
over a length of the order of the superconducting coherence length, = fivp/A, 
and therefore the electron-hole conversion occurs over this distance. Thus, the phase 
shift during Andreev reflection is sensitive to the profile A(x) of the superconducting 
order parameter over a distance of the order of a coherence length near the interface. 
Calculation of this profile involves self-consistent treatment of the order parameter: 
however, the properties of the surface states can be qualitatively understood by 
replacing the exact form of A(a;) by a step function A(a;) = AO(a;). 

2.2. Surface states 

Let us now consider the formation of superconducting bound states in an INS quantum 
well with a specular vacuum/metal (IN) interface surface. We will first consider a 
two-dimensional s-wave superconductor with anisotropic order parameter which has 
a perfectly transparent normal metal/superconductor (NS) interface with a normal 
layer of width L, as shown in figure 1. According to quasiclassical arguments, the 
bound state corresponds to a closed quasiparticle trajectory, and the bound state 
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energy is given by the Bohr-Sommerfeld quantization condition which requires the 

total phase accumulated during one cycle to be equal to an integer number of 2tt. 
The quasiparticle trajectory in figure 1 (corresponding to electron motion in the +ky- 
direction) consists of an electron segment which includes a single reflection at the IN 
interface, and of a hole segment which backwards retraces the electron trajectory. The 
total accumulated phase consists of two parts: (i) a contribution from the Andreev 
reflections: —7 — x \n the ^-direction and —7 + X in the ^-direction {9 is the angle 
between the electron trajectory and the surface normal and = tc — 9 \s the angle 
after specular normal reflection at the surface) , and (ii) a contribution from the phase 
/3 accumulated during propagation through the normal region, /3 = 2L{k^ — k^) + /?o, 
where the first term corresponds to the ballistic motion and the second to the reflection 
at the NI interface. The electron/hole wave vectors are hk'^''^ = 2m{Ep cos^ 9± E). 
The Bohr-Sommerfeld quantization condition can then be written down as 

-{7 + l)TSx + PiE) = 2nn, (1) 

where Sx = X~X^ and where T^x corresponds to quasiparticle trajectories in the ±ky 
directions. 

In the case of a uniform and isotropic order parameter (7 = 7; X — x)j 
equation (1) reduces to 

-27 + /3 = 2n7r (2) 

and leads to the spectral equation 

i; = ±Acos^^. (3) 

According to this equation, surface states always exist, even at a simple IS interface 
{L = 0) due to the normal reflection phase shift /3o 7^ in the superconducting 
electrode (Wendin and Shumeiko, 1996a). However, in the presence of an INS well, 
L ^ 0, the ballistic contribution to /3, 2L{k'^—k^) « ALE/Uvf cos 9, usually dominates, 
and the dispersion equation (3) takes the form found by de Gennes and Saint-James 
(1963), 

E ^ E 2L 

— = ±cos — -. (4) 

A A^ocos6' ^ ' 

Generally, the electron-hole dephasing factor (3{E) is an odd function of the 
energy and turns to zero at = 0. This is a result of the fact that the difference 
between electrons and holes vanishes at the Fermi surface. Due to this property there 
is no E' = solution of equation (3), and the density of surface states in s-wave 
superconductors is zero at the Fermi level; see figure 2(a). Moreover, due to the 
presence of time reversal symmetry in the problem, the surface levels always appear 
in pairs symmetrically positioned around the Fermi level. 

The number of surface states for each quasiparticle trajectory is determined by 
the length of the quasiparticle path, which is given by the normal layer thickness and 
the propagation angle 9: for L < tt^Qi equation (4) implies that for normal incidence 
{9 = 0) there is a single pair of surface states close to the energy gap edge. A new 
pair of surface states appear at the gap edge when l/cos0 = ttS^q/L. Increasing the 
trajectory angle further pushes these states down into the gap, and a new pair of gap 
edge states appears every time 1/ cos 9 is increased by tt^^q/L. Integration over the 
angle 9 yields the surface density of states plotted in figure 2(a), cf. reference (de 
Gennes and Saint- James, 1963). 
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In anisotropic s-wave superconductors, different gap values before (|A|) and after 
(|A|) normal reflection at the surface (7 7^ 7) will modify the phase shifts of the 
Andreev reflections and change the energy of the surface state. Therefore the surface 
states depend on how the surface is oriented relative to the crystal axes. The angle 
integration will then to a large extent wash out the peaks in the density of states (DOS) 
found in figure 2(a). Roughness of the interface will lead to further suppression of the 
peaks in the DOS and to the opening of a gap in the spectrum (Melsen et al., 1996). 
In this review we will however not discuss rough interfaces in any detail. It is also 
worth mentioning that for anisotropic superconductors, the total DOS within the 
energy interval Amin < E < Amax, contains, in addition to the bound surface states, 
contributions from the continuum states. 

The surface states are two-fold degenerate, the two states corresponding to 
different signs of the electron momentum (±fcy) parallel to the surface. In the presence 
of a magnetic field parallel to the surface, this degeneracy is removed due to the 
Meissner current: the spectrum calculated above holds in the frame of reference of the 
condensate. In the laboratory reference frame the levels are shifted by the Doppler 
term p^vp sinO, E ^ E + pgVp sinO, where ps is the superfluid momentum of the 
Meissner current. This shift has opposite signs for electrons moving along and against 
the Meissner current. 

2.3. d-wave surface states 

In the d-wave superconductors, the order parameter is angle dependent. 



where a describes the orientation of the order parameter relative to the surface normal; 
see figure 3. Consequently, the gap (order parameter) is different before (A) and 
after (A) normal reflection. As long as A and A have the same sign, i.e. the same 
phase, X = X; the situation is not qualitatively different from the anisotropic s-wave 
case. However, the situation dramatically changes if A and A have opposite signs, 
as illustrated in figure 3: in this case, the phases x said x = X + of these order 
parameters differ by (5x = X ^ X = ^'^i which leads to a modification of equation (2), 



This equation always has the root E = since ^{E = 0) = ^{E = 0) = 7r/2 
and (}{E = 0) = 0. This means that for the corresponding angle of incidence, a 
surface state will appear at zero-energy: this is the so-called midgap state, MGS 
(Hu, 1994). Clearly, this MGS is degenerate with respect to ±ky. It is easy to see 
that for the a = ±7r/4 orientations, an MGS exists for every quasiparticle trajectory. 
— 7r/2 < 9 < 7r/2, while for the a = 0,n/2 orientations there are no MGS at all. For 
arbitrary a, MGS exist within the window 7r/2 — a < |^| < 7r/2 -|- a. The formation 
of MGS is a robust property of d-wave superconductors: it requires only different 
signs of A and A and does not depend on the possible difference in their absolute 
values. Moreover, the formation of MGS is not sensitive to the spatial variation of 
A(,'r) (see Appendix A. 4). This insensitivity is an important property, because in 
d-wave superconductors the gap is suppressed by the surface for most orientations of 
the crystallographic axes. 

There is work in the literature addressing the questions of how surface roughness 
(Golubov and Kupriyanov, 1999), impurities (Poenicke et al., 1999), and nanofaceting 
of the surfaces/interfaces (Fogelstrom et al., 1997; Tanuma et al., 1998) affect the MGS. 



A{e) = Aocos[2(0-a)], 



(5) 



-{j + j)+TT + l3{E) = 2mT. 



(6) 
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However, in this review we will focus on clean ballistic structures with specularly 
reflecting interfaces, and we refer the reader to the original literature (e.g. the 
abovementioned papers and references therein). 

The major effect of taking gap suppression into account will be the presence 
of additional bound states in the vicinity of the gap edge. Fully self-consistent 
calculations have been performed by e.g. (Nagato and Nagai, 1995; Buchholtz 
et al., 1995a; Buchholtz et al., 19956; Matsumoto and Shiba, 1995a; Matsumoto 
and Shiba, 19956; Matsumoto and Shiba, 1996; Barash et al., 1997). Qualitatively, 
the effect of the surface levels at finite energy can be understood by considering a 
rectangular SNI well of size L ~ ^o- For the sake of illustration, let us consider the 
a = 7r/4 orientation for which an MGS exists for every trajectory angle (see figure 3). 
We then have |A| = |A|, and equation (6) reduces to 



The resulting density of bound states for a normal metal length L = 1.5^o if' shown 
in figure 2(b). Note that the density of surface states in this figure is only the part of 
the density of states with \E\ < A^axi similar to the anisotropic s-wave case; there is 
also a considerable contribution from the extended states of the bulk due to the gap 
nodes. 

2.4- Impurity induced MGS 

The above discussion was focused on MGS at planar surfaces. However, it is clear 

that MGS must exist at any defect which produces quasiparticle scattering between 
states and corresponding to different signs of the order parameter: cracks, twin 
boundaries, impurities, etc. Thus, a direct measurement of the DOS at zero energy 
by scanning tunneling microscopy (STM) in the vicinity of an impurity can test the 
d-wave symmetry of the order parameter. Such experiments were recently done by 
Yazdani et al. (1999), Hudson et al. (1999) and Pan et al. (2000a). We shall discuss 
these experiments later, at the end of Sec IV, which is devoted to a discussion of 
current transport through MGS. 

2.5. Paramagnetic effect of MGS 

Another possibility to detect the impurity induced MGS was demonstrated in the 
experiment by Walter et al. (1998) (figure 4), where the response to an external 
magnetic field was studied. In this experiment, the MGS were introduced in a c- 
oriented YBCO film by creating multiple quasi-lD in-plane defects (ion tracks) by 
ion bombardment. By varying the orientation of the ion beam, and therefore the ion 
tracks, with respect to the order parameter, it was possible to control the density of 
MGS. In the experiment, the temperature dependence of the magnetic penetration 
depth A(T) was measured for different orientations of the defects (figure 4). It was 
found that for (110) orientation A(T) has a minimum below 15K and increased with 
decreasing temperature, while for (100) orientation the low-temperature upturn was 
much smaller. 

This effect can be explained within the above picture by studying the 
paramagnetic response of MGS to an applied magnetic field (Higashitani, 1997). As 
was mentioned in section 2.2, the surface states undergo splitting in the presence of 
a magnetic field. The energy of an individual MGS becomes Emgs{^) = PpVsy sin 6, 




|A| ^ocos^' 



E 2L 



(7) 
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where Vgy is the superfluid velocity along the surface of the defect. Since the MGS 
arc localized over a distance much smaller than the penetration depth (^q ^)-. 
the nonhomogeneity of the superflow can be neglected. The MGS current per ab- 
plane along the surface is jmgs = ^VFSva.6\'il)MGs{x)\^nF[EMGs{())]- The total 
current of all MGS at the defect is calculated by integrating over x and also over 
ky, —kp < ky < kp- At zero temperature, T = 0, only the negative energy levels 
{ky < 0), which carry paramagnetic current, are occupied, which leads to a total 
current Jmos — -^eEp /'^ttH at the defect [for (110) orientation]. By dividing this 
equation by d, the average distance between the defects, we get the average density 
of the paramagnetic MGS current. At finite temperature, the positive energy levels 
carrying diamagnctic current also become occupied, which leads to a reduction of 
the paramagnetic efi'ect. When T ^ PfVs, expansion of the Fermi function and 
integration over ky yield a 1/T-dependence for the average MGS current density, 
juGS = —{eEp /QTT'h,d){ppVs/T). This paramagnetic current is to be added to the 
paramagnetic current of the quasiparticles in the bulk of d-wave superconductor, 
jn = —eVsUniT), where HniT) is the density of quasiparticles which linearly depends 
on the temperature. nn{T) ^ {ksT/ A)n, at ksT <C A (n is the electron density). 

This interplay of the bulk quasiparticle current and the MGS current produces 
the low-temperature anomaly of the magnetic penetration depth observed in the 
experiment (Walter et al., 1998). The effect is most pronounced for (110) orientation, 
and disappears for (100) orientation of the defects because of the absence of MGS in 
the latter case. In the experiment, the effect did actually not disappear for the (100) 
orientation, but was a minimum. The finite value was explained to be due to MGS at 
internal interfaces in the film, e.g. grain boundaries. 

2.6. Splitting of MGS 

Although the d-wave symmetry of the gap necessarily leads to the formation of 
MGS, the resulting large density of states at the Fermi level is actually energetically 
unfavorable: any mechanism able to split the MGS and produce a gap in the spectrum 
will lower the energy of the system (Sigrist, 1998; Fogelstrom and Yip, 1998). Different 
mechanisms splitting the MGS have been proposed in the literature: (1) splitting due 
to a subdominant component of the order parameter which forms a local complex d±is 
gap parameter at the surface, (2) self-induced Doppler shifts, and (3) formation of spin- 
density wave near the surface. We shall here discuss the first two mechanisms in greater 
detail, while for the third scenario, we refer to the original literature (Honerkamp and 
Sigrist, 1999; Honerkamp et al., 2000). 

2.6.1. d± is order parameter For orientations where the d-wave gap is suppressed 
and at low temperatures, a subdominant gap of s or dxy symmetry may appear near 
the surface if a weak secondary interaction (additional to the one responsible for 
the d-wave superconductivity) is present (Matsumoto and Shiba, 1995 a; Matsumoto 
and Shiba, 1995&; Kashiwaya et al., 1995; Matsumoto and Shiba, 1996; Fogelstrom 
et al., 1997; Tamuiia et al., 1999; Zhu et al., 1999). In a self-consistent picture, the 
complex combination A^; ^ ^ ± iAg (or Ad .-^ .-^ ± iA^^ ) is energetically favorable. 
The MGS is sensitive to the presence of this subdominant gap; to demonstrate this, 
let us consider the order parameter A^ .-^ .-^ + iAg with the d-wave component having 
a = tt/4 orientation. We now define a relative phase angle Xrei = arctan(As/A(i). 
In this case, the scattering phase shift during Andreev reflection has the additional 
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contribution +Xrel foi" fc^ > 0, and —Xrel for ky <Q. In a symmetric case with 7 = 7, 
we get the following modification of equation (6) at low energies: 



The effect of this ±fcj^-dependent contribution from the relative phase is to shift the 
MGS away from the Fermi level to the edges of the subdominant gap, 



and the ±.ky degeneracy of the MGS is lifted. This splitting of the MGS is a result 
of the broken time-reversal symmetry in the + iAs state. The MGS splitting 
drives the instability (lowers the energy of the system) which appears at a critical 
temperature determined by the strength of the subdominant interaction. Each bound 
state carries current parallel to the surface in the direction given by ky, which means 
that at low temperatures, when only the state below the Fermi energy is occupied 
[the one labeled by +ky for > in equation (9)], a net surface current will be 
present. This spontaneous current manifests the broken time- reversal symmetry of 
the order parameter (Sigrist, 1998). The TRSB state is two-fold degenerate: the gap 
parameter can be cither d + is 01 d — is, corresponding to the two possible directions 
of the spontaneous surface current. 

The spectrum of surface states in d + is superconductors has a close similarity 
to the spectrum in a B-phasc p-wavc superconductor. In the light of the above 
discussion, it is therefore straightforward to understand the results obtained by 
Buchholtz and Zwicknagl (1981). The B-phase gap parameter in a two-dimensional 
p-wave superconductor can be written A^^ = A((t cos6 + i sin 9) (cr = ±1 for down/up 
spins), see e.g. the book by VoUhardt and Wolfle (1990). This gap function is 
similar to the d + is case above, but with a propagation-angle-dependent magnitude 
of the imaginary component zAsin^. Thus, there are surface bound states also 
for p-wave superconductors; however, the bound state is located at the Fermi level 
only for quasiparticles propagating perpendicular to the surface {9 = 0), since the 
gap is real (and changes sign at the normal reflection event) only for this particular 
propagation angle. Integrating over the injection angle, we find that the density of 
states is finite for all subgap energies, and there is no zero energy peak (Buchholtz 
and Zwicknagl, 1981). We conclude that the density of boimd surface states of p-wave 
and d-wave superconductors are quite different, although the mechanism of the bound 
state formation is the same. 

2.6.2. Doppler shifts In connection with the discussion of the paramagnetic effect, 
Higashitani (1997) discovered that the MGS are unstable, even in the absence of 
any subdominant interaction (pure d-wavc order parameter), due to a possibility of 
creation of spontaneous magnetic fields. As we discussed in the previous section, a 
splitting of MGS yields a surface current at low temperature. This current creates a 
magnetic field which is screened by a supercurrent counterflow which, in turn, further 
splits the MGS due to the Doppler shift effect discussed in section 2.2. The effect 
saturates when the magnetic energy compensates for the energy of the Doppler shifted 
MGS. Thus, a time-reversal symmetry breaking (TRSB) state spontaneously occurs, 
with a surface current due to the split MGS and a self-consistent screening current 
upholding the Doppler shifts at a phase transition temperature ksTTRSB/ Aq = 
(l/6)(6/A) (Honerkamp and Sigrist, 1999; Barash et al., 2000; Lofwander et al., 2000) 
which is much smaller than the critical temperature of the type-II superconductors 



27 =F 2xrei +TT = 2mr . 



(8) 



Emgs = -(sgnfcy) A, 



(9) 



(6 « A). 
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3. dc Josephson effect 

In this Section we shall discuss how, from a theoretical point of view, MGS influence 
the dc Josephson effect. First, we develop the concept of Andreev bound states and 
explain how they in many cases determine the properties of Josephson junctions. 
Then we discuss some particular junction orientations which illustrate some effects 
peculiar for Josephson junctions of rf-wave superconductors: large Josephson currents 
proportional to ^/D (D is the jimction transparency), 1/T dependence (T is the 
temperature) of the Josephson critical current, 7r-junction formation, ^ tt junction 
crossovers with temperature, and spontaneous time-reversal symmetry breaking 
(TRSB). 

3.1. Andreev states: coupled surface states 

When two superconducting electrodes are connected to form a junction, the 
superconducting surface states form coherent bound states of the entire junction, 
Andreev bound states. There is an obvious analogy with a double well system. The 
structure of the Andreev state is illustrated in figure 5. Like the surface states, 
the Andreev state consists of a combination of electron and hole wave functions 
with opposite electron momenta kx (momentum ky parallel to the specular surface 
is conserved). The analogy with the double well system allows us to expect that 
the Andreev levels will be split and shifted with respect to the surface levels, the 
shift depending (i) on the transparency of the tunnel barrier separating the electrodes 
and (ii) on the superconducting phase difference across the junction. Taking into 
account the phase dispersion of the energy {dE/d(j)) of the Andreev levels, one may 
talk about Andreev bands E{(j)) with width (dispersion) proportional to the junction 
transparency D. However, if the surface states at the two sides of the junction have 
equal energies, the coupling becomes resonant. In this resonant case, the splitting 
of the levels, and therefore the width of the Andreev band, will be particularly large 
and proportional to VD (Wendin and Shumeiko, 1996a). As we will see, the width of 
the Andreev band determines the critical Josephson current, and the Andreev band 
dispersion determines the Josephson current-phase relation. One might expect the 
surface states to be degenerate in junctions fabricated with the same superconducting 
material; however, even in this case the degeneracy can be lifted due to several types 
of asymmetries, e.g. misorientation of the order parameters, difference in spatial 
modulation of the the order parameters, etc. 

The derivation of the spectral equation for the Andreev states is more involved 
than the derivation of the spectrum of surface states. The properties of the Andreev 
levels can therefore not be strictly deduced from simple quasiclassical arguments but 
rather require solving the Bogoliubov-de Gcnnes (BdG) equation (dc Gcnnes. 1989). 
However, the spectral equation can be presented in terms of the earlier introduced 
phase shifts -jl (7/?) and /3l {0r) in the left (L) and right (R) electrodes of the 
junction, and the superconducting phase difference 4> = xl — XR- We here choose to 
define the Andreev reflection phase 7 for negative A as 7 = arccos(ii^/|A|) + tt. The 
spectral equation then becomes 

cosa = i?(6') cos 6 + £)(6') cose, (10) 

where 2a = -71, - Ti, - 7ij - Ta + /3l + I3r\ 26 = -71, - + 1r + 1r + I3l - Pr; 
2c = -71, + 71, - 7ij + 7fl - 2(f). 
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Equation (10) describes hybridization of the surface states. This is obvious 
if one lets the transparency D{6) 0, in which case one recovers equation (6) 
for decoupled surface states. However, equation (10) can also be interpreted as 
hybridization of Andreev states present in the corresponding perfectly transparent 
junction {D = 1) due to normal electron back scattering at the insulating barrier. 
The spectral equation for Andreev states in perfectly transparent junctions actually 
foUows from quasiclassical arguments: 

-IL -1R + pL+ pR - 4> = 2mr, 

-IL -iR + h+PR-<P = 2n7r, (11) 

where the first and second lines correspond to the decoupled 6 and 6 segments in 
figure 5. 



3.2. Josephson current 

Since Andreev reflection is accompanied by charge transfer through the NIS interface, 
the Andreev states are able to carry current. The current through the surface states 
is blocked in the direction perpendicular to the surface due to complete quasiparticle 
reflection from the insulator. In the NIS junction, this insulating interface becomes 
transmissive, and the current may flow through the Andreev states. There is a general 
relation between the current through the Andreev state and the phase dispersion of 
the energy of the Andreev state, 

This equation can be derived directly from the BdG equation (Shumeiko et al., 1997) or 
deduced from the thermodynamical equation I = {2e/fi)dF/d(p by using a microscopic 

expression for the jimction free energy (Beenakker, 1991). The total Josephson current 
of the bound states per unit surface area, per a6-plane, flowing perpendicular to the 
surface is 

(13) 



(•••) = / 



where the angle brackets denote integration over the trajectory angle 9, 

fTT/2 

decose. (14) 

rr/2 

The Josephson current in d-wave junctions strongly depends on the orientation 
of the order parameters in the electrodes with respect to the junction interface. 
The whole picture is quite complex, however one can distinguish three qualitatively 
different cases: (i) no MGS at any side of the junction, (ii) MGS at both sides of 
the junction, and (iii) MGS at one side of the junction. For arbitrary orientation all 
three cases may occur for different sectors of the trajectory angle 0, and the total 
current is a weighted sum of currents of the three types. In a pure form, the first case 
is realized in do/c^o junctions where the situation is similar to s-wave junctions: no 
suppression of the order parameter, and the two Andreev bound states which dominate 
the Josephson current are situated close to the gap edges and have dispersion of order 
D. The second case is realized in a pure form in d^ j^j dj^^ junctions. In this case, 
the central phenomenon is the resonant coupling of MGS which produces large current 
proportional to ^/D with rather special properties. We choose to illustrate the third 
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case by the do /c?7r/4 junction. The pecuUarity of this case is the time-reversal symmetry 
breaking, TRSB. 

In many works, e.g. (Sigrist and Rice, 1995), it has been assumed that tunneling 
occurs only for normal incidence, ^ = 0. This is motivated by the discriminating 
effect of the tunneling cone (Wolf, 1985). However, in order to include the effects of 
MGS it is very important actually to take a wide tunneling cone D(9) into account, 
since tunneling through MGS only occurs at finite angles of incidence. Ignoring this 
point, one might conclude that the Joscphson current is zero, or at least small, for the 
d^/4/d^/4 orientation since A(0 = 0) = 0. However, this will be true only for a cone 
very strongly peaked around ^ = 0, since the MGS contribution for finite cone angle is 
resonant. Clearly, the exact functional dependence of the current on orientation will 
be sensitive to the actual shape of the cone and not only to the value of the junction 
transparency {D) . Throughout this review we use a J-function model for the barrier. 
In this case the cone is quite wide, D{0) ^ cos^ 6 [sec e.g. (Bruder, 1990)], and the 
effects of MGS are somewhat overemphasized. An exception is the treatment of the 
do/d-^/i orientation below, where we instead follow reference (Lofwander et al., 2000) 
and use a wide barrier model which has D{6) ^ expi—KO"^), where k determines the 
width of the cone and is a function of the barrier height and width (Wolf, 1985). 
However, even in this case the MGS contribution will dominate at low temperature, 
as discussed below. 

3.2.1. d^/i/dT^/i orientation: 0-junction with large Ic The order parameter for this 
orientation is strongly suppressed by the surface; however, the qualitative properties 
of MGS in junctions with low transparency D ^ 1 are not sensitive to the detailed 
profile of the order parameter, which allows us to consider the rectangular well model. 

For this high-symmetry orientation, the four gaps in equation (10) have equal 
magnitudes, |Aj:,| = = |Afl| = \Ait\. The tt shift at both sides is picked up at 
9 (Andreev reflecting against a negative lobe of the order parameter). The spectral 
equation takes the form 

cos {2'y-f3) = -R{e) + D{e) cos {(}>), (15) 

which, for small energies <C A and negligible dephasing (/3 = 0), has the solution 
(Riedel and Bagwell, 1998) 

£;± = ±|A|/Dcos^, (16) 

plotted in figure 6(a). The current of these two MGS is calculated via equation (13). 
At zero temperature, only the level below zero energy is populated while the level 
above is empty, and the current is 

JMGS = ^ (|A|\/iJ^ sin^ sgn ^cos0 . (17) 

The critical (maximum) Josephson current is proportional to Vd (Tanaka and 
Kashiwaya, 1996a; Riedel and Bagwell, 1998) and is much larger than the Ambegaokar- 
Baratoff critical current in non-resonant s-wave tunnel junctions, which is proportional 
to D (Ambegokar and Baratoff, 1963). In figure 6(b), we plot the current-phase 
relation at low temperature. Note that the \/^-behavior is a general phenomenon 
connected with resonance coupling (Wcndin and Shumciko, 1996a). 

To understand the basic temperature dependence of the current, we notice that 
the two MGS energies E± carry currents in opposite directions, j+ = —j-, and 
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therefore with increasing population of the upper level at large temperature, the 
current carried by MGS will decrease, 

2ekF /dE+ ,^ , dE_ ,^ , 

2ekF / dE+ , E+ \ 
= r-^ ( tanh — !- ) . (18) 

The scale of the temperature dependence of the current is set by the well-defined energy 
difference of the two-level MGS system. When Ej^ -C /cbT, the expansion tanh a; w x 
leads to a 1/T dependence of the current (Tanaka and Kashiwaya, 1996 a; Barash 
et al., 1996; Riedel and BagweU, 1998), 

In tunnel junctions, the current carried by MGS dominates at low temperature 
fcsT < because the continuum state current is proportional to the junction 
transparency D. Due to suppression of the gap at the interface, there are also Andreev 
states near the gap edges which may carry large currents due to degeneracy of the 
corresponding surface states. However, the not cm-rent of these states is also of order 
D since they form pairs far below the Fermi level which are almost equally populated. 
The currents carried by the two bound states of the pair flow in opposite directions 
and cancel each other. Thus we conclude that the critical current Ic scales with 
temperature as 1/T for low temperatures, as clearly seen in figure 7 (the solid line). 
In the limit T — > 0, the /c(r) oc l/T dependence approaches a maximum value 
7c(0) a \i\\\fD, c.f. equation (17). 



3.2.2. d^/4/c?-7r/4 orientation: ir-junction with large Ic From a technical point of 
view, this case is very similar to the previous one, the difference being that the tt shift 
at the right side is now picked up at 9 rather than at 6 (because the positive and 
negative lobes of the order parameter have changed places). This results in a phase 
shift (p ^ <p + n in equation (15) which now reads 

cos (27 -f3) = -R{e) - D{e) cos (^) , (20) 

with the MGS solution (Riedel and Bagwell, 1998) 

i;± = ±|A|\/^sin^. (21) 

The MGS band again has a large dispersion proportional to Aa/D; however, the phase 
dependence is tt shifted compared to the rf,r/4/c^7r/4 junctions [sec figure 6(c)]. The 
current-phase relation will therefore also be 7r-shifted [c.f. figure 6(d)] which puts 
the minimum of the junction energy at a phase difference equal to tt instead of 0, 
corresponding to a 7r-junction. This fact was used in the tri-crystal ring experiments, 
although the theoretical idea behind the experiment put forward by Geshkcnbcin et al. 
(1987) and later on by Sigrist and Rice (1995) was based on symmetry arguments and 
not on the microscopic picture outlined here. 

Again, the critical current has the low-temperature anomaly, Ic{T) 1/T and 
approaches a maximum value /c(0) oc |A|-\/!d at zero temperature, (dashed line in 
figure 7). To emphasize the 7r-junction behavior, we let the critical current be negative 
in this case. 
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3.2.3. da/d-a orientation: to n -junction cross- over We are now able to discuss 
general properties of the Joseplison current in symmetric junctions d^/d^a] sec 
(Barash et al., 1996; Tanaka and Kashiwaya, 1996a; Tanaka and Kashiwaya, 1997; 
Tanaka and Kashiwaya, 1999) for a detailed discussion. In these junctions, MGS exist 
only within certain limited intervals of the angle while for other angles, the Andreev 
level spectrum is similar to the one in s-wave junctions and the corresponding currents 
follow the Ambegaokar-Baratoff law. At low temperature, the MGS will dominate and 
the equilibrimn phase difference is, as in the d-^j^jd^T^j^ case, at (6 = tt. However, 
at larger temperatures this contribution will be suppressed and the non-MGS current 
may dominate. This s-wave type current will give a minimum of the junction energy 
at (/) = 0. Hence, with increase of the temperature, the junction will undergo a 
transition from a 7r-junction to an ordinary junction (Barash et al., 1996; Tanaka and 
Kashiwaya, 1996a). The critical current is zero at the transition temperature where 
there is detailed balance between the MGS current and the s-wave like current. We 
illustrate this behavior with the d^/g/d-Tr/s junction, the dot-dashed line in figure 7. 

3.2. J^. d^j/^ -\- is / d_.,^ii + is orientation: junction re-entrance Development of a 
subdominant s-wave order parameter near the junction interface leads to an interesting 
modification of the low-temperature anomaly of the critical current. As we earlier 
discussed, in the presence of a complex order parameter -I- lAg the MGS are split 
into two states Emgs = where the sign depends on the sign of the ky wave 

vector. Thus the Andreev spectrum presented in figure 6 will move either up or down 
(for opposite signs of ky) as soon as the subdominant gap appears and grows with 
decreasing temperature. Let us consider the limit As(0) » |A|-\/^ where the effect 
is most pronounced. In this case, the Andreev bands do not cross the Fermi level at 
T = 0; the lower pair of levels is then fully occupied, while the upper one is empty 
and therefore does not contribute to the current. This leads to an uncompensated 
t/-current along the interface. Large x-currents carried by the occupied levels flow in 
opposite directions and cancel each other, leading to a small (for D <C 1) residual 
current proportional to the first power of the junction transparency D. This residual 
current has positive sign, similar to the common situation in s-wave SIS and SNS 
junctions (Wendin and Shumeiko, 1996a; Samuelsson et al., 2000), and therefore the 
TT-j unction property is lost. Thus the Josephson current undergoes a crossover from 
TT-junction to 0-junction when the temperature decreases, and the critical current 
has a non-monotonic temperature dependence and turns to zero at the crossovcir 
point (Tanaka and Kashiwaya, 1998), see figure 8. Experimental observation of such 
reentrant temperature dependence of the critical current would be a direct evidence 
of the subdominant order parameter. 

3.2.5. do/dT^ii orientation: -periodicity For the ^0/1^^/4 orientation, there are MGS 
only at one side of the junction corresponding to the non-resonant case. The spatial 
shapes of the order parameters are different at the two sides of the junction which 
also leads to non-resonant coupling of the non-MGS surface states. To understand 
the type of bound states present in this junction, we solve the spectral equation for 
the injection angle 9 = n/S when all four d-wave gaps involved have equal magnitudes 
|A(±7r/8)| = A0/V2. The spectral equation then takes the form [D = £)(7r/8), A = 
A(7r/8)] 

sin(27 -I3) = ±D sin {(/)) (22) 
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where the different signs correspond to quasiparticle states with ^ky. The solution 
for /? = is given by 

£2 = ^ (^1 ± Y^l-£)2sin2<^^ . (23) 

There are two Andreev bands for a given ky, the MGS band and the band at the gap 
edge (Riedel and Bagwell, 1998). The dispersion of the MGS band is proportional to 
D at low transparency, 

Emgs = -sgnky^Dsm(t> + 0{D^), (24) 

while the dispersion of the Andreev band at the gap edge is much smaller, proportional 
to D^. Equation (24) is a good approximation for all injection angles in the low 
transparency limit when the MGS are well inside the gap, close to zero. 

For the orientations discussed in the previous sections, the MGS spectrum is 
degenerate with respect to ±ky: in other words, there is a pair of MGS, E±, for every 
ky. In the present case, however, this degeneracy is lifted as emphasized in figure 9(a). 
Inspection of the MGS wave functions for phase differences between and tt shows 
that, at large distance from the interface, the E- wave function decays only for +ky, 
while for —ky it grows exponentially. For the E+ wave function the situation is the 
opposite: the wave function decays for —ky and grows exponentially for +ky. This 
difference between the two time-reversed states labeled by ±ky is a consequence of 
the sign difference between the d-wave lobes and is the key point in understanding the 
possibility of spontaneous time- reversal symmetry breaking (see below). 

The current in the c?o/'^7r/4 junction is 7r-periodic, as opposed to the usual 2tt- 
periodicity, as was first theoretically predicted in (Yip et al., 1990; Yip, 1993; Tanaka 
and Kashiwaya, 1997; Zagoskin, 1997; Ostlund, 1998) and later on measured in 
experiments by Ilichev et al. (1999). For low temperature the MGS band dominates 
the current-phase relation plotted in figure 9(b) and the bound states in equation (23) 
give Ic{T) cx D^IAp/fc^T. However, for large temperatures Ic will be negligibly 
small already for temperatures of the order |A|£), in contrast to the d^^ j d^^ j case 
where this happens for temperatures of the order |A|\/^. This difference is a result 
of the larger energy scale given by the large dispersion of the bound states (cx VZ?) 
in the resonant d^/i/d^n case compared to the much smaller dispersion (cx D) in the 
non-resonant d^/dT^n case. 

3.2.6. do/d^^ii orientation: TRSB The current through the junction is zero for the 
phase differences </> = n7r/2. However, the minimum junction energy is achieved for 
0eq = ±7r/2, since the MGS energy is the lowest for these phase differences. Due to 
the uneven occupation of MGS with opposite signs of ky, the current parallel to the 
interface is finite even in equilibrium, see figure 9(c). For example at (j)eq = 7r/2, only 
MGS with positive ky are occupied and a surface current is flowing in the positive 
y-direction. The opposite happens for = — 7r/2. Remembering that the state with 
—ky is obtained from the +ky state by the time-reversal operation, and noting that 
there; is no way of taking the time-reversed state to the original one by adding an integer 
number of 27r to (j), we must conclude that time-reversal symmetry is spontaneously 
broken at the do/'^7r/4 junction. This fact is manifested by the finite surface current, 
which produces a magnetic field at the junction. 

The Josephson energy gain for = ±7r/2 is partly compensated by the 
energy cost of setting up the Meissner screening currents which ensure that the 
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spontaneous magnetic field vanishes in the bulk superconductors far from the junction. 
By calculating the thermodynamic potential for the junction system (Lofwander 
et al., 2000) it can be shown that the TRSB state really is favorable. If the orientation 
is not exactly do/d^/i, the equilibrium phase difference (j)eq is shifted from ±7r/2. 
However, the TRSB effect (i.e. 4)cq 7^ or tt) remains for deviations up to about 10° 
(Yip, 1997; Lofwander et al., 2000) from the do/d^i^ orientation. 

Since there are screening currents flowing, we may expect Doppler shift effects. 
For transparencies D ^ Co/A, the dispersion of the MGS with phase difference, 
given by equation (24), is much larger than the Doppler shifts; these can therefore 
be neglected and the picture outlined above holds. However, for D <C Co/A the 
Doppler shifts dominate, and the problem is formally the same as MGS at a free 
surface discussed in section 2.6. 

The above discussion of the rfo/rf,r/4 junction, including the crossover to the 
doppler shift driven TRSB surface instability, closely follows the one presented 
in the paper by Lofwander et al. (2000). The possibility of TRSB in purely d- 
wave junctions was first pointed out by Yip (1995) for the pin hole geometry, 
which is equivalent to having D = \. Many people, using different theoretical 
tools, have considered the do/dTr/i junction: by using Ginzburg-Landau theories 
(Sigrist et al., 1995; Ostlund, 1998) or quasiclassical Green's function techniques 
(Yip, 1995; Yip, 1997; Fogelstrom et al., 1998; Fogelstrom and Yip, 1998), or by solving 
the Bogoliubov-de Gennes equation (Huck et al., 1997; Zagoskin, 1997; Zagoskin and 
Oshikawa, 1998; Lofwander et al., 2000). However, the final results are qualitatively 
the same: 7r-periodic c;m-rc;nt-phase relation and the possibility of spontaneous time- 
reversal symmetry breaking. There is however one important difference: the absence 
(as above) or presence of a subdominant component of the order parameter at the 
junction. We conclude that broken time-reversal symmetry at a cio/c^7r/4 Josephson 
junction will not give any evidence for the presence of a subdominant component, 
since the d-wave order parameter alone enforces broken symmetry. 

In experiments reported by Kirtley et al. (1995a) and Mannhart et al. (1996a) 
fractional fluxes were found at grain-boundaries with the do/c?OT/4 orientation. If time- 
reversal symmetry is broken as described above, spontaneous surface currents offer a 
possible explanation for these experiments, see (Sigrist, 1998). However, Mannhart 
et al. (1996a; 19966) proposed an explanation in terms of randomly distributed facets 
along the grain-boundary, where each facet may have equilibrium phase difference 
or TT depending on whether zero, one or two lobes with negative sign of the rf-wavc 
order parameter point towards the junction. When the phase difference varies along 
the grain-boundary, formation of fractional fluxes may be energetically favorable. It 
should be noted that this model assumes a simple sin dependence of the current- 
phase relation, and does not take into account the formation of MGS at the junction. 

3.2.7. do /do junctions with resonant levels in the barrier There is an interesting 
similarity between the Josephson effect in the d^/4/cJ^/4 junctions and in do /do 
junctions with resonant states in the tunnel barrier. In the presence in the tunnel 
transparency of a narrow Breit-Wigner resonance close to the Fermi level, the 
Andreev level spectrum has the form found in the s-wave case (Beenakker and van 
Houten, 1991; Wendin and Shumeiko, 19966), 




(25) 
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where Eq is the position and V is the width of the resonance, -BcF <C A. If the 
resonance is exactly at the Fermi level, Eq = Q, the spectrum becomes 



which is the spectrum of an SNS-junction {D = 1) with | A| = F, having dispersion F 
instead of |A|. At zero temperature, the current is then given by 



The temperature dependence of the critical current in this case is then also 
IciT) oc 1/T for UbT ^ F. Thus, while interpreting experimental results it is useful to 
have in mind that the low-temperature anomaly of the critical current is a fingerprint 
of a resonant state at the Fermi level rather than indication of the d-wavc symmetry, 
although this resonance (MGS) could be due to the d-wave symmetry of the order 
parameter. 

3.3. Critical Josephson current in grain boundary and ramp-edge junctions 

An important class of experimental studies involves the determination of the critical 
Josephson current in current-biased artificial grain boundary (GB) junctions da^ /dan 
by varying the orientation angles and of the order parameters (Gross, 1994; 
Delin and Kleinsasser, 1996; Ivanov et al., 1998; Aric ot al., 2000). 

Consider a gedanken experiment where we may freely rotate the order parameters 
of the left (ai) and right {an) superconducting electrodes, and also independently 
vary the transparency D of the interface barrier (I) (figure 5), creating ideal rfat/rfas 
junctions: how will the Josephson current develop in different situations? 

Starting with an ideal do /do GB junction, there is no MGS, and the Josephson 
current oc £> is due to non-resonant transport through the interface barrier. If we now 
rotate the order parameter symmetrically away from a = 0toQ; = 7r/4, keeping D 
fixed, the da/d-a the junction will develop MGS, and the critical current will then 
cross over from D to y/D dependence because of resonant transport through MGS. 
The resulting current enhancement can become very large (" giant Josephson current" ) 
if the transparency is low, <C 1, and the IcRn product will be enhanced by 1/VD 
over the Ambegaokar-Baratoff limit characterizing do /do tunnel junctions. 

In reality one must rotate the order parameter by rotating the lattice. This 
necessarily leads to an interface barrier, the strength of which increases with 
misorientation angle (up to 2a = w/A for tetragonal symmetry). Assuming that 
the thickness of the barrier increases with misorientation angle, the result will be 
exponentially decreasing transparency and narrowing of the timneling cone. With 
decreasing transparency, the low-temperature region where MGS are important will 
decrease {ksT oc a/ {D)Ao for resonant MGS). These effects may seriously reduce the 
weight of the MGS contribution. Therefore, in real experiments, the role of MGS will 
depend on the detailed properties of the junction. 

It should be noted that if we continue to rotate the structure with the order 
parameters to reach the d^/4/d_^/4^ limit, there is no longer any lattice; mismatch, 
and therefore no interface barrier. The full presence of MGS therefore cannot be 
achieved in this type of experiment: to really observe the effects of MGS, one must 
fabricate d^/4^/ 1 /d_^/4 junctions with separately controlled tunnel barriers. 

In real high-Tc GB junctions, the IcRn product decreases strongly with increasing 
misorientation angle. In the absence of MGS, one would expect this product to 



E = ±Tcos{(t)/2), 



(26) 




(27) 
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be constant, but it is experimentally often found that /c-Rjv oc 1/pN, where is 
the normal resistance per unit area (Gross, 1994). An explanation, due to Gross 
(1994), is that Coulomb blockade prevents resonant Josephson current but allows 
resonant normal current transport through the levels in the barrier. We suggest 
that this behavior can also be explained in terms of a constant spectral density of 
normal levels with resonance widths F inside the barrier. From the discussion in 
section 3.2.7, we know that for a resonance with energy Eq <^ F, equation (27) gives 
a resonantly enhanced Josephson current / oc F, while for resonances with i?o S> F 
the current is negligibly small. The contribution to the total current then comes from 
a spectral region of width F around E = 0, with the results / oc F^. The normal 
current, on the other hand, is resonantly enhanced in the entire spectral region, giving 
IcRn oc F oc 1/pn- 

In a recent paper, Arie et al. (2000) have investigated the properties of the 

Josephson current in ramp-cdgc junctions with a/j — = 7r/4 overall misoricntation 
angle, varying the orientation of the interface (and therefore an, q;l) to vary the role 
of MGS. In particular, the dependence of the critical current on temperature has 
been experimentally investigated and fitted to d-wavc theoretical results for short d/d 
junctions with high transparency, D=0.5-0.8, which implies that the dispersion of the 
MGS is a large fraction of the gap. To get agreement with experiment for T — > 0, Arie 
et al. (2000) introduced life-time broadening of the MGS (describing effects of rough 
interfaces and inelastic scattering) via an imaginary part of the energy, which strongly 
reduces Ic{T) and the IcRn product at low temperatures, changing the sign of the 
curvature and leading to saturation of Ic{T). The final good qualitative agreement 
between theory and experiment is consistent with the presence of life-time broadened 
MGS in a high-transparency junction. 

It should be noted, however, that even if the transparency for superconductive 
transport seems high, an estimate of the experimental average transparency, based 
on the junction resistance (0.2 — 0.8f2) and area (number of planes and transport 
channels), gives a fairly low value, around D = 10~^ — 10"'*. We suggest that 
the simplest explanation for this difference is that transport takes place through 
highly transparent point contacts in the junction barrier. In this picture, the 
temperature dependence of the critical Josephson current becomes a measure of the 
true transparency (£) ps 1), which together with the normal resistance also becomes a 
measure of the effective transport area of the junction. 

For a more definitive analysis it is desirable to compare experimental and 
theoretical results for a larger set of misorientation angles and transparencies. To 
this end we have calculated Ic{T) for high transparency, (D) = 0.68 [Z = 0.5, 
as in the original fit by Arie et al. (2000)], neglecting MGS broadening. We have 
considered the high-symmetry orientations do /do and 'i7r/4/<i-7r/4) as well for the 
intermediate orientation d_ 150/^300 considered by Arie et al. (2000). The result is 
shown in figure 10. In the case of do /do orientation (no MGS) we find that Ic{T) has a 
weak negative curvature, while in the case of dj^/4/d_j^/4 orientation (MGS is present) 
we find a weak positive curvature of Ic{T)- Thus, within our model the positive 
curvature is characteristic for the presence of MGS, while the negative curvature is 
characteristic for the absence of MGS. The strong positive curvature in the case of 
the rf_i5o/rf3oo orientation in figure 10 may look a bit surprising: in fact, only at low 
temperatures docs the l/T-dcpendencc signify MGS; at higher temperatures it results 
from a competition between the nearly-linear MGS contribution and contributions of 
opposite sign from the continuum and gap-edge states. Thus, theoretically, positive 
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curvature of 7c {T) at low temperatures is a fingerprint of MGS. However, such a picture 
needs to be verified by systematic experimental and theoretical investigations, and the 
picture may be complicated. For example, the recent studies of GB junctions by Tafuri 
et al. (1999) show clear positive Ic{T) curvatures for nominally do/do orientations, and 
such behavior is also evident in the review by Delin and Kleinsasser (1996), perhaps 
due to orientational disorder like meandering and faceting. 

Arie et al. (2000) also studied the magnetic field dependence of the critical 
Josephson current, finding evidence for sin(2(/)) behavior. This might be evidence 
for do/c?OT/4 symmetry, but may also be due to the high transparency (saw- tooth- like 
dependence of the /(^) current-phase relation). 

We conclude that so far, the role of well-defined MGS is difficult to identify with 
certainty in studies of the critical Josephson current in d/d grain-boundary or ramp- 
edge junctions, although Arie et al. (2000) most likely have taken an important step in 
the right direction. In order directly to observe the presence of MGS one may have to 
use low-transparency tunnel junctions and utilize energy-resolving probes, e.g looking 
for MGS in the form of characteristic features in I-V characteristics of voltage biased 
junctions. 

4. Voltage-biased N/d junctions 

Since the discovery of tunneling in NS junctions (Giaever, 1960), and the explanation of 
the phenomenon on the basis of the tunnel Hamiltonian model (Bardeen, 1961; Cohen 
et al., 1962), it is commonly accepted that the tunneling conductance essentially 
measures the superconductor density of states (DOS). Hence, tunnel spectroscopy 
provides a direct means of detecting the MGS, which should show up as a large 
conductance peak centered at zero bias. This zero-bias conductance peak (ZBCP) 
has indeed been observed; see (Geerk et al., 1988; Kashiwaya et al., 1995; Covington 
et al., 1997; Ekin et al., 1997; Alff et al., 1997; Wei et al., 1998; ApriH et al., 1998; Alff 
et al., 1998a; Alff et al., 19986; Sinha and Ng, 1998; Suzuki et al., 1999; Krupkc and 
Deutscher, 1999; Aprih et al., 1999; Wang et al., 1999; Covington and Greene, 2000) 
and references therein. See also the recent reviews by Alff and Gross (1999) and 
Kashiwaya and Tanaka (1999; 2000). 

We start this section by discussing the theory of tunneling through MGS. We will 
discuss two different transport mechanisms: direct tunneling, and tunneling involving 
Andreev reflection. The tunnel model is usually exploited for calculation of direct 
tunneling, while Andreev reflection is usually described by means of scattering theory. 
We will compare the two approaches and show that depending on the relation between 
the junction transparency and the MGS intrinsic broadening, the MGS resonance may 
appear either in the single particle current or in the pair current. In the end we discuss 
tunnel experiments on NS junctions and also the recent STM investigations of impurity 
bound states. 

Theoretical model 

In the tunnel model, the angle-resolved local density of states (LDOS) pL and pR 
are calculated on the left(L) and right(R) sides of the junction, and then the tunnel 

formula 




(28) 
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is employed to calculate the current. The total current is obtained by angle integration 
I{y) = (/(y,^)), and the conductance is then found by differentiation with respect 
to the bias voltage V. Since the density of states (DOS) in the normal metal is 
approximately energy independent near the Fermi surface, in the zero temperature 
limit the conductance will be reduced to {Dps{eV)), where ps is the LDOS in the 
superconductor . 

For the a = orientation of the order parameter, there is no MGS and the 
angle resolved LDOS is the same as in the bulk ps oc \E\/ ^/W^\^Q{\E\ - |A|) 
[see e.g. references (Matsumoto and Shiba, 1995c; Barash et al., 1995)]. The tunnel 
formula gives a conductance which in the zero temperature limit is simply the bulk 
DOS corrected by the angle dependent transparency D = D{6) of the tunnel barrier, 

G{V)ozId , 9(ey-|A|)). (29) 

For other orientations, the conductance will measure LDOS at the surface rather 
than the bulk DOS. For the a = 7r/4 orientation, neglecting gap suppression and 
surface roughness, the LDOS is given by (Matsumoto and Shiba, 1995c) 

p,{e, E) (X Q(^\E\ - |A|) + TTimE). (30) 

The high spectral weight at the gap edges has moved down to zero energy and 
formed the MGS, and the tunnel conductance will contain a peak at zero voltage. 
The delta-function singularity is usually removed by introducing broadening via an 
imaginary part t] of the energy representing some relaxation mechanism. This damping 
T] determines the width of the MGS resonance. 

This interpretation of the conductance spectra in terms of single-particle 
(quasiparticle) current, see figure 11(a), assumes that all states in the electrodes form 
reservoirs, i.e. their relax;ation times are assumed short enough to ensure equilibrium. 
This assumption is reasonable for the c^xtcnded states of the continuum, while for the 
MGS, lying deep within the superconducting gap, it requires special consideration. If 
the MGS are true bound states without coupling to the bulk of the superconductor, 
i.e. intrinsically very narrow, then current transport through the MGS is only possible 
via the Andreev process, see figure 11(b). The equation for the current then reads 

I{V)=h{V)+l2{V), (31) 

where /i and I2 are the single-particle (quasiparticle) and two-particle (Andreev) 
currents respectively. This case of narrow MGS is most easily studied within the 
scattering approach (Landauer, 1957; Biittiker, 1986; Imry, 1986). 

Blonder et al. (1982) were the first to use scattering theory to calculate the 
conductance of an ,s-wavc NS junction. Later on, Bruder (1990), Kashiwaya et al. 
(1995) and Xu et al. (1996) generalized the theory to the d-w&ve case; see also (Tanaka 
and Kashiwaya, 1995; Kashiwaya and Tanaka, 1996; Tanaka and Kashiwaya, 19966). 
At zero temperature, in the low-transparency limit, we have for the a = 7r/4 
orientation 

h{e,E) = D '^^'-^' e{\E\ - |A|), (32) 

and the two-particle current takes the form 

2D2 
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The total current is obtained by integration over angles and energies 

I{V) (X 



J dE[nF{E)-nF{E + eV)][h{e,E)+l2{0,E)]^ (34) 



The two-particle current contains a Breit-Wigncr resonance of width AD/2\/R, the 
MGS, and I2 oc AD as soon as eV > AD. Hence, the two-particle current is resonantly 
enhanced near zero bias (it is of the same order as the single-particle current), which 
produces a zero-bias conductance peak. The current-voltage and conductance-voltage 
relations calculated within the scattering approach are plotted in figure 12. 

Consider now the ratio q between the life time due to leakage to the normal 
reservoir through the barrier, T5 cx Ti/AqD, and the life time due to inelastic relaxation, 
Tr oc h/r], which couple the MGS with the bulk superconductor. 

Clearly this life-time ratio will determine whether the MGS will participate in current 
transport via the single-particle process (for large damping; q ^ 1), or via the two- 
particle process (for small damping; q <^1). 

The fingerprint of Andreev transport is the so-called excess current defined as 
(Artemenko et al., 1979) 

/exc= lim [IsiV)- iNiV)], (36) 

where Is and In are the currents in the superconducting and normal states 
respectively. Following the work of Blonder et al. (1982), it became common practice to 
estimate the efficiency of Andreev rcfiection from the magnitude of the excess current. 
Since in the present case I2 ~ D, one would expect that also I^xc ~ D. However, it 
can be shown that the single-particle current gives a large negative contribution to 
the excess current: this cancels the enhanced two-particle current and leaves an excess 
current which for all orientations is proportional to in the low-transparency limit, 
just like in the low-temperature superconductors. This implies that the excess current 
will not be a true measure of the importance of the Andreev reflection process as soon 
as there are resonances (like the MGS) in the junction. 



4-2. Experiments 

Experimentally, the zero-bias conductance peak (ZBCP) has been studied for many 
years: it was observed by Geerk et al. already in 1988. In the early years, when 
the rf-wave symmetry of the order parameter had not been established, the peak was 
explained in terms of a spin-fiip scattering model developed by Appelbaum (1966) 
and Anderson (1966). Following the paper by Hu (1994), attention was drawn to 
the MGS scenario. One of the first good fits between experiments and the drwave 
theory (MGS, figure 13) was presented by Kashiwaya et al. (1995). Today, a number 
of experiments have given quite convincing evidence that the MGS picture offers the 
best explanation. Here we shall discuss some experimental findings regarding the 
influence of (1) magnetic fleld and (2) orientation and disorder on the MGS and 
ZBCP. Yet another class of experiments have been performed recently, namely STM 
measurements of MGS bound to point impurities. These experiments will be discussed 
in section 4.2.3 below. 
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4.2.1. Magnetic field The Doppler shift effect mentioned in section 2.2 was confirmed 
experimentally by Covington ct al. (1997) and may be regarded as the first successful 
attempt to seriously discriminate between the MGS and the spin fiip scattering 
scenarios for the ZBCP. The effect was further studied by Krupke and Deutscher 
(1999) and ApriH ct al. (1999), see figure 14. The ZBCP was observed to spHt linearly 
with H up to Hci, where vortices enter the superconductor and the shift saturates 
(Fogelstrom et al., 1997). It was pointed out that the spin-flip model could be ruled 
out since the Landc g-factor involved in the Zeeman type splitting in this case had to 
be very large in order to fit the experiment. 

In the zero-field limit i? — > 0, splitting of the ZBCP was still observed by 
Covington ct al. (1997) and Krupke and Deutscher (1999) and also by Kashiwaya 
et al. (19986) and Lesueur et al. (1999). This was interpreted as a sign of spontaneous 
time-reversal symmetry breaking at the surface and evidence for a local presence of a 
subdominant s-wave component, c.f. section 2.2 above. However, a spontaneous (zero- 
field) splitting of the ZBCP is not always seen (Aprih et al., 1999; Yeh et al., 2000). 
Moreover, there are other possible explanations for the splitting than the one involving 
a subdominant ,s-wavc component of the order parameter; see the discussion of split 
MGS in section 2.6. Also, if for any reason particle-hole symmetry is broken in the 
superconductor, the MGS may also be split. For the identification of the mechanism 
responsible for the spontaneous splitting of the ZBCP the same arguments as for 
the identification of the origin of the ZBCP are relevant: one needs to measure the 
response to some controllable perturbation; e.g. the magnetic field response of the 
ZBCP discussed above confirmed that MGS arc responsible for the ZBCP. Since such 
experiments are lacking for the spontaneously split ZBCP, we feel that at the present 
time the existence of a TRSB d + is surface state remains controversial. 

4.2.2. Orientation As explained above, if we probe an a = 7r/4 surface (i.e. [110] 
surface) a large ZBCP should be present, while if we probe an a = surface (i.e. [100] 
or [010] surfaces) the ZBCP should be absent. In some well controlled measurements 
(Alff et al., 1997; Wei et al., 1998; Wang et al., 1999; Iguchi et al., 2000) it has been 
verified that the ZBCP has this orientation dependence (see figure 15). However, 
because of junction inhomogeneities like faceting (Fogelstrom et al., 1997), the local 
orientation at the junction may not be well defined. For example the planar junction 
experiments performed by Covington et al. (1997) showed equal weight of the ZBCP 
for all orientations of the electrodes. 

Aprili et al. (1998) studied broadening effects on the ZBCP due to impurities 
introduced into the junction via ion bombardment. The result suggests that disorder 
introduces dccohcrcncc which removes the MGS resonance. Aprili et al. (1998) argue 
that if the spin-flip mechanism were responsible for the ZBCP, the disorder should not 
have such a large impact, indicating that the MGS are responsible for the ZBCP. 

Figure 15 shows that the peak height decreases as 1/T while the width seems to 
remain roughly constant, independent of temperature. This information is obtained 
by assuming that the peak is positioned on a background in the form of a temperature- 
dc;pendcnt gap. In that case, the FWHM is roughly the distance between the two "fix 
points" in figure 15. This behavior implies that the density of MGS (the area of the 
peak) decreases with increasing temperature, which can not be explained by thermal 
smearing. The question of temperature-independent width has been addressed by 
Walker and Pairor (1999) in terms of the effect of rough interfaces and surface umklapp 
scattering. 
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4-2.3. Impurity induced MGS Bound states induced by impurities have been found 
experimentally by different groups (Yazdani ct al., 1999; Hudson ct al., 1999; Pan 
et al., 2000a). In the most recent experiment by Pan et al. (2000a) the LDOS 
was measured by STM in the vicinity of a single impurity (Zn) replacing a Cu 
atom in the superconducting plane of BSCCO (figure 16). It was found that the 
LDOS is enhanced at low energy and it has an anisotropic form in real space. 
This observation is in qualitatively good agreement with theoretical predictions; see 
(Balatsky ct al., 1995; Balatsky and Salkola, 1996; Salkola ct al., 1996; Flattc and 
Byers, 1998; Tsuchiura et al., 1999; Tsuchiura et al., 2000; Zhu et al., 2000a; Zhu et 
al., 20006; Haas and Maki, 2000) and references therein. It was predicted that a point- 
like impurity in a d-wavc superconductor is able to form a resonant state within the 
superconducting gap. If global particle-hole symmetry is present, then in the unitary 
limit (strong impurity potential), the energy of the resonance approaches the Fermi 
level, and the resonance width tends to zero, the MGS becoming a true bound state. 

We can understand this effect by considering a defect with a hard specular 
boundary the size of which is much larger than the coherence length, as discussed by 
Matsumoto et al. (1997). In this case, the defect surface can be locally approximated 
by a planar surface, and our arguments in favor of the presence of the MGS holds. 
Since different parts of the defect boundary have different orientations with respect 
to the crystallographic axes and continuously change from a = to a = 27r, there 
are different amounts of trajectories containing local MGS. Consider for example a 
circular defect, as illustrated in figure 17. Then for the (110) and (-110) directions 
the local angle is a = 7r/4(l + 2n) which corresponds to a maximum amount of MGS 
trajectories, while for (100) and (010) orientations, the local angle is a = n7r/2 which 
corresponds to a minimum amount of MGS trajectories. Hence, the LDOS at low 
energy (Matsumoto ct al., 1997) is cross shaped in space with tails along the nodes. 
The density of states decreases far from the defect because the MGS are localized near 
the scatterer over a distance of the order of the coherence length. 

In fact, the above arguments rely entirely on quasiclassical dynamics of 
superconducting quasiparticles, and can therefore be extended (see Appendix A. 4) 
to defects smaller than the coherence length but larger than the Fermi wave length 
(Adagidcli ct al., 1999). However, in high-Tc materials this difference is not that 
important since the coherence length is rather small. 

At the time of writing, an issue under debate is whether the impurity MGS are 
localized or extended, sec e.g. (Zhu ct al., 2000c). In a large crystal there are many 
impurities, or inhomogeneities, and hence many MGS, which in principle can form 
an impurity band similar to impurity bands in semiconductors. In simplified terms, 
the MGS wavefunction can leak out along the gap nodes (the cross-shaped DOS 
described above); if wavefunction overlaps between neighboring impurity states are 
significant, an impurity band may form. Also experimentally, the issue is controversial. 
Thermal conductivity measurements (Taillefer et al., 1997; Behnia ct al., 1999) on 
Y123 and BSCCO indicate the existence of extended quasiparticle states, while other 
investigations (Hussey et al., 2000) on under-doped Y124 indicate the opposite. 

5. Voltage-biased d/d junctions 

In this Section we study the ac Josephson effect in superconducting junctions. We 
begin by outlining the microscopic theory of charge transport in voltage biased 
Josephson junctions. Then we discuss how the MGS resonance appears in this case. In 
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particular, we show that depending on the hfetime ratio q introduced in the previous 
section, the MGS resonance appears either in the pair current at the gap voltage 
= A (for small intrinsic broadening, g <C 1), or in the single particle current both 
near zero voltage and at the gap voltage (for large intrinsic broadening, 1). In 
the end we discuss some experiments. 

5.1. Theoretical formulation 

Generalization of the previous theory of transport through N/d junctions to d/d 
junctions [beyond the tunnel model; equation (28)] is non-trivial due to the Josephson 
effect. In the presence of a bias voltage, the phase difference over the Josephson 
junction becomes time-dependent, dtcj) = ^eV/h, (Josephson, 1962), and a quasiparticle 
incident on the junction will see the junction as a scatterer whose properties 
periodically change in time with the Josephson frequency loj = 2eV/h. In the limit 
of small applied voltage, eV A, the time-dependent current can be understood via 
adiabatic arguments, in terms of a time-dependent phase, letting (j) (j){t) ~ 2eVt/h 
in the equation for the dc Josephson current (Averin and Bardas, 1995; Bratus' 
et al., 1997). In the general case, a full quantum mechanical consideration is necessary 
for inelastic scattering by the time-dependent scatterer, taking into account transitions 
between the incident quasiparticle energy E and an infinite set of sideband energies 
En = E + neV (n are integers) . The time-dependent current can then be decomposed 
into a spectrum of frequency components 

oo 

I{V,t)= J2 Im{V)e^"'^^\ I.m{V) = lUV), (37) 

m=— oo 

where the m = component is the time averaged current and the m 7^ components 
correspond to the ac oscillations with basic frequency ujj = 2eV/h. Interestingly, the 
very early observation by Esteve et al. (1987) of the ac Josephson effect was taken as 
one of the first proofs for the existence of Cooper pairs in the high-Tc superconductors. 

Within the adiabatic picture for the d-wave junctions, the 7r-periodic current- 
phase relation in the rfo/rf7r/4 junction (see figure 9) will result in an ac Josephson 
oscillation with twice the usual frequency, 2ujj. This was predicted by Yip et al. 
(1990; 1993) and Zagoskin (1997) using the adiabaticity arguments, and it was 
recently confirmed by a full calculation of equation (37) by Lofwander et al. (1998). 
Possibly this effect was observed in experiments by Early et al. (1993), although the 
authors presented an alternative explanation in terms of parallel junctions in the grain 
boundary. 

The time-independent (m = 0) current, from now on called the dc current, has 
contributions from all sidebands, — 00 < En < 00, and at low temperature it has the 
form 

00 

Iac{V) = Y.C{V). (38) 

Remarkably, this expansion for the dc current has two different physical 
interpretations, which in fact are equivalent, namely multiparticle tunneling (MPT) 
and multiple Andreev reflection (MAR). A consistent theory unifying the concepts 
of MPT and MAR has been presented by Bratus' et al. (1995) and Shumeiko et al. 
(1997). 
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The first interpretation, multiparticle tunneling (MPT), is particularly useful 

in the tunnel limit, D 1. It was originally suggested by Schrieffcr and Wilkins 
(1963) using perturbative tunnel model calculations, and it was improved by Cuevas 
et al. (1996). According to this interpretation the n = 1 current, if'^iV), is due 
to singlc-particlc tunneling processes, and it corresponds to the quasiparticle current 
in a standard tunnel model calculation. The n = 2 current, /|'^(y), is due to two- 
particle tunneling processes which have smaller probability, ~ D^, and so on. In 
s-wavc junctions, energy conservation prohibits the n-particle tunneling process for 
neV < 2 A. Thus, each n-particle current has a threshold at voltage eVn = 2A/n and 
the magnitude of this current is proportional to £>". However, resonances due to the 
presence of surface states or resonances in the barrier may lower the exponent and 
enhance the intensity of the process (Johansson et al., 1999a; Johansson et al., 19996). 

The second interpretation, more useful for transparent junctions, is based on 
the scattering theory and involves multiple Andreev reflections (MAR) (Klapwijk 
et al., 1982; Arnold, 1987; Bratus' et al., 1995; Averin and Bardas, 1995). According 
to this interpretation, the current Ii''{V) is due to direct scattering (acceleration) 
across the energy gap, see figure 18(a), while the current /|'^(y) is due to the Andreev 
reflection processes, in the same way as in the NS junction, equation (31), depicted 
in figure 18(b). For smaller voltage, the Andreev reflected hole may be Andreev 
reflected once more when reaching the left superconductor, this time as an electron 
[see Fig. 18(c)], giving the current I^'^{V), and so on. The multiple Andreev reflections 
form, in principle, an infinite series of currents, equation (38). The magnitude of 
each current is proportional to _D" (in the absence of resonances), since each process 
contains n coherent passages through the barrier. 

The equivalence of these two interpretations — MPT and MAR — is due to the 
fact that at every Andreev reflection one pair of electrons is transmitted through the 
SN interface, and therefore the MAR process with n Andreev reflections transfers a 
charge (n + l)e, corresponding to the n + 1-particle current (Arnold, 1987; Shumeiko 
et al, 1997). 

In s- wave junctions, the onsets of the n-particle currents at eV = 2A/n, together 
with peaks at the same voltages in the higher order currents due to the DOS peaks 
near the gap edges, is the explanation for the so-called subharmonic gap structure 
(SGS), see (Ludoph et al., 2000) and references therein. 

In d-wave superconductors, the sub-gap structure (Hurd, 1997; Barash and 
Svidzinsky, 1997; Hurd et al., 1999; Lofwandcr, Johansson, Shumeiko, Wendin and 
Hurd, 1999; Samanta and Datta, 1998; Yoshida et al., 1999; Lofwander, Johansson 
and Wendin, 1999) is very different for two reasons: (1) the multiparticle current 
thresholds are washed out due to the gap nodes, (2) the presence of MGS changes the 
resonances in energy space. For the do /do orientation (no MGS in the junction), the 
single-particle current dominates the I — V characteristics because of the gap nodes. 
Higher order processes (n > 2) are very small corrections and can be neglected; thus, 
sub-harmonic gap structure at eVn = 2A/n is suppressed in the I-V characteristics 
in rf-wave junctions, see Fig. 19(a). When MGS are present in the junction, new 
resonances appear. In particular, the two-particle process becomes resonant (see 
below) and must be taken into account. Although the MGS may enhance also higher 
order processes than n = 2, these resonances are always very weak compared to the 
two-particle process. For these reasons we will in the following only discuss the single 
particle and two particle currents. 
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We start our discussion of the I — V characteristics of junctions with MGS by 

presenting in figure 19(b)-(d) results of full numerical calculations of equation (38) 
for junctions where the right electrode has d-wave symmetry with orientation an = 
IT /A and the left counter electrode has s-wave symmetry (s/fi^/4 junction), d-wave 
symmetry with = orientation {do/d-^/^ junction), and d-wave symmetry with 
OiL = 7r/4 orientation {d^i^/d^ii junction), respectively. In the last junction 
case we introduced a small imaginary part (to model broadening) to the energy 
in order to reveal the current peak near zero bias. If the MGS are strongly 
coupled to the reservoirs (large intrinsic broadening, q ^ 1), the current- voltage 
characteristics in figure 19 are dominated by single particle tunneling, and they can 
be qualitatively understood by applying tunnel model arguments. Indeed, when the 
right superconductor is a d-wave superconductor with the 7r/4 orientation, the local 
density of states is given by equation (30); the MGS delta function 5{E) then produces 
a term (Samanta and Datta, 1998) 

lt%V) (X {D{e) [np{-eV) - np^ ppiO, -eV)) . (39) 

Within this model, we understand figures 19(b)-(d) as simply the s-wave DOS, the 
angle integrated d-wave DOS, and finally the midgap peak. Note that the current itself 
rather than the conductance is proportional to the LDOS pp of the left superconductor. 
For voltages eV > A terms not included in equation (39) become important and for 
eV ^ A the current is proportional to the applied voltage. 

When the MGS are intrinsically very narrow (q <C 1), the resonance will not 
appear in the single-particle current, instead the resonance appears in the two-particle 
current (compare the discussion of the A''/d junction in section 4). We exemplify this 
with the .5/^71-/4 and d^/4/d^/4 junction cases by plotting the single- and two-particle 
currents separately along with the total currents in figure 20. In the s/(i^/4 junction 
case, the I-V relations (total current) for the g » 1 and 5 <C 1 cases are similar, 
although in the g <C 1 limit the peak at the s-wave gap voltage is due to the two- 
particle current; the resonant process is shown in the inset. This situation is equivalent 
to the one in the N/d case discussed in section 4. However, in the d^^j^jd^j^ junctions, 
which have MGS at both sides of the junction barrier, the situation is different. The 
process of MGS-to-MGS tunneling is absent in the small q limit. Thus there is no 
current peak at small voltage in this case [compare figure 19(d) and figure 20(d)]; 
instead the resonance shows up in the two-particle current at the gap- voltage eV = A, 
see figure 20(f); the resonant process is shown in the inset. Thus, for this orientation 
it should be possible experimentally to distinguish the two cases, large or small q, by 
observing where the resonance appears: either near zero voltage {q ^ 1) or at the gap 
voltage (g <C 1), respectively. 

5.2. Comments on experimental results 

A large number of current-voltage characteristics (IVC) measurements have been 
performed with grain boundary and ramp-edge junctions (Delin and Kleinsasser, 
1996). In this Review we will not make any serious attempt to analyze these 
experiments beyond what is stated in the original papers. Instead we will make some 
comments on a few specific issues. 

Measurements of IVC of d/d junctions are difiicult to analyze. Besides 
complications with the samples (see the discussion of measurements of the critical 
Josephson current in Section 3.3), the presence of the Josephson effect makes the 
sub-gap voltage region hard to interprete. 
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Because of the Josephson current, an essential issue concerns how to measure 
IVC. In order to measure quasi-particlc current and conductivity, and to observe the 
MGS zero-bias peaks, one must work with voltage bias, which by definition avoids the 
problem with the dc Josephson current. However, this is hard to do experimentally, 
because of the low junction impedance; in practice, current bias is therefore usually 
employed. 

Alff et al. (1998&; 1998a) have recently reported IVC measurements on symmetric 
da/d-a HTS grain boimdary jimctions with large misoricntation angles, in the range 
30 — 45°. These junctions show very small or vanishing Josephson current, and the 
conductance dl/dV vs V shows a zero-energy structure quite similar to NIS junctions 
(see figure 15), with similar temperature and magnetic field dependences. Alff ct 
al. (19986; 1998a) conclude that they in fact observe NIS-like behavior, although 
they interpret the zero-energy structure in the dl/dV vs V characteristic as a ZBCP 
resulting from an SIS convolution of density of states with strongly broadened MGS. 

The similarity of the ZBCPs of the NIS and SIS junctions is striking, and we would 
like to consider the following tentative explanation. For an SIS junction to behave like 
an NIS junction, it is necessary that there are very strong relaxation effects in the 
barrier region, the barrier acting as a reservoir, which results in decoupling of the two 
SI and IS interfaces. The SI/IS junctions will then be in series and each will be biased 
by eV/2; as a result, there will be a zero-bias onset of the current (without peak) as in 
the NIS junction, but the maximum gap should occur at 2A. The average gap voltage 
Vg could be well below this value. Whether this agrees with the experimental results 
remains an open question since there is no independent measurement of the gap in 
the experiment. 

Alff et al. (19986; 1998a) rather consider their results to signify that the junction 

is in the tunneling regime, SIS, albeit tunneling is considered to proceed through 
resonant states in the barrier. Let us then assume that we can describe the junction 
as a "tunnel" junction S"r'S, where the "F'-region cannot be regarded a reservoir. The 
Josephson current could still be suppressed for several reasons: the junction could be 
long in comparison with the coherence length, or suppressed by inelastic interactions, 
or there could be fluctuations and noise causing the two sides of the junction to 
become incoherent. In a symmetric situation with narrow MGS on both sides of the 
junction, as described before, the MGS DOS peaks will then be convoluted, leading 
to a current peak, i.e. to a peak in the IVC and not in the conductance. In this view, 
the conductance zero-energy anomaly in the da/d^a junction describes the voltage 
derivative of the MGS DOS, and should not be directly compared with the MGS DOS 
of the NIS junction. 

Note that if the MGS is very broad, say T ^ 0.3Ao, the MGS peak in the 
calculated IVC will merge into the rising background, and the result may appear 
similar to the NIS case; nevertheless, in our opinion, one should still look for the MGS 
DOS in the IVC, and not in the conductance. However, with such extreme broadening 
the MGS may act as normal reservoirs in the gap, and the whole situation may have 
to be reconsidered. 

In conclusion, we believe that further progress depends on close interaction 
between experiment and theory to critically examine models and results in order to 
produce internally consistent pictures and new critical experimental tests. 
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6. Other experiments 

6.1. c-axis tunneling 

For a long time, the presence of a Josephson current in experiments on c-axis tunnel 
junctions between YBCO and Pb has been regarded as rather mysterious, because if 
the order parameter of YBCO has pure rf-wavc symmetry, the first order Josephson 
coupling should vanish, only leaving a small residual second order current (see e.g. 
references (Sun et al., 1994; Kleiner et al., 1996; Lesueur et al., 1997). This conclusion 
however relics on the assiimption that the c-axis current injection, from e.g. an STM 
tip, results in a spherically symmetric ab-plane injection current. Only under such 
circumstances is the overlap between the s-wave and the rf-wave order parameters zero. 
Considering spherical injection, the orthorhombic structure of YBCO indicates that 
there may be, in addition to the d-wave order parameter, a small s-wave component 
present which could be responsible for the first order Josephson coupling observed in 
the experiments. However, twinning spoils this explanation since neighboring twins 
should have different signs of the s-wave component, and the average Josephson 
current should be close to zero. Kouznetsov et al. (1997) studied tunneling along 
the c-axis from Pb into YBCO containing one single twin boundary. It was shown 
that indeed a d± s order parameter is realized in YBCO and that the sign between 
the rf-wave and s-wave components is opposite on the two opposite sides of the twin 
boundary. However, also in experiments on junctions between Pb and highly twinned 
YBCO (Kleiner et al., 1996), a first order Josephson couphng was found. Furthermore, 
in experiments on c-axis BSCCO-Pb junctions (Mossle and Kleiner, 1999) a finite first- 
order Josephson coupling was observed and interpreted in terms of a small s-wavc 
component also in BSCCO. Recently an interesting possibility was brought forward 
by Rae (2000), who pointed out that Pb is a rather anisotropic s-wave superconductor 
and showed how this asymmetry on the s-wave superconductor side can explain at 
least the experiments on BSCCO. To test this idea, a thorough investigation of the 
Josephson coupling between isotropic s-wave superconductors and BSCCO should be 
performed. 

6.2. Electron- doped cuprates 

Very recently, Tsuei and Kirtley (20006) found the half- flux quantum also in tricrystal 

rings of the electron-doped cuprates NCCO and PCCO. This provides very strong 
evidence for d-wave superconductivity in both these materials. This conclusion 
is supported by the power-law dependence of the penetration depth found by 
Prozorov et al. (2000). However, in transport experiments (Ekin et al., 1997; Alff 
et al., 19986; Kashiwaya et al., 1998a; Kleefisch et al., 1998; Alff et al., 1999) ZBCPs 
were never found. Instead, more or less V-shaped spectra were obtained. This 
absence of ZBPC was taken as evidence for an anisotropic s-wavc order parameter in 
NCCO, but now turns out to be in disagreement with the very recent tricrystal ring 
experiments (Tsuei and Kirtley, 20006). Therefore, at this point, a reexamination of 
the transport properties of the electron-doped materials seems necessary in order to 
check the presence of MGS also in these materials. We would like to point out that 
the experimental results obtained by Alff can be understood within the MGS scenario 
if the MGS in these materials are intrinsically very narrow [small q, equation (35)]. 
However, also in the normal-metal-superconductor (N/d) experiments on NCCO by 
Kashiwaya et al. (1998a) the ZBCP was absent, which can not be explained by a 
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narrow MGS only. 

7. Summary and concluding remarks 

Midgap states (MGS) are consequences of the angle-dependent sign change of the order 
parameter and provide direct evidence for d-wave symmetry. Transport measurements 
which detect the MGS resonance at the Fermi energy for particular orientations 
of the plane of the junction with respect to the crystallographic axes will provide 
crucial evidence for rf-wave superconductivity. The MGS phenomenon is part of a 
general problem of superconducting surface and interface states and their relation 
to the Andreev states, which determine the transport properties of superconducting 
junctions. 

Surface midgap states (MGS) in d-wave superconductors are robust features 
which result in a large number of observable effects: large Joscphson current, low- 
temperature anomaly of the critical current, 7r-junction behavior. ^ tt junction 
crossover with temperature, zero-bias conductance peaks, paramagnetic currents, 
time reversal symmetry breaking and spontaneous interface currents, and resonance 
features in subgap currents. 

Some of these consequences of MGS have been verified experimentally (zero- 
bias c;on(luctancc peaks, paramagnetic currents and impurity-induced states), while 
some remain controversial (time-reversal symmetry breaking and spontaneous 
surface/ interface currents), or so far unobserved (large Josephson current and the 
1/T low-tcmpcrature anomaly of the critical current in low-transparency jimctions). 

A number of problems need to be clarified, e.g. the character and strength of 
relaxation processes and the width of the MGS. In particular, does absence of MGS 
imply absence of rf-wave (and presence of s-wave), or does it imply d-wave and strong 
decoherence? 

The predicted low-temperature anomaly associated with MGS, 7c oc 1/T, remains 
to be confirmed experimentally. The lack of experimental evidence may be due to 
faceting of the grain-boundary (Hilgenkamp et al., 1996): different parts of the junction 
may have local orientations supporting large MGS currents, but with opposite signs 
for different facets. Consequently, the low-temperature anomaly associated with the 
MGS may be averaged out for a macroscopic junction. The obvious remedy is to work 
with junctions with interfaces with well-defined orientation. 

The true nature of HTS grain-boundary (GB) and step-edge Josephson junctions 
remains, in our opinion, a mystery. It is difficult to get a consistent picture. In recent 
work on low-impedance step-edge junctions, a reasonable description of the Josephson 
current could be obtained assuming high transparency, while at the same time the 
low junction resistivity suggested that the transport is through high- transparency 
spots, point contacts or pinholes, in the interface barrier. In GB junctions at 
moderate misoricntation angles a similar description may be appropriate. At large 
misorientation angles, however, the GB junctions has been characterized as an SIS 
tunneling junction with resonant states in the barrier. However, it is difficult to judge 
whether the experimental results imply coherent tunneling or incoherent transport 
through the junction. Here coevolution of experiment and theory seems to be necessary 
to obtain deeper understanding. 

Other predictions which need to be verified experimentally include negative 
differential conductance in sub-gap structure (needs a voltage-bias setup), and 
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spontaneous magnetization in time-reversal symmetry breaking (TRSB) states at 
do/c?7r/4 junctions. 

In a number of interesting recent experiments, a picture has emerged that 
the superconducting order parameter of the cuprates might not have pure dx2_y2 
symmetry but rather have a complex order parameter with d,j.2_y2 + id^y symmetry, 
where the dxy component is small compared to the dx2_y2 component. STM studies 
of vortex cores in BSCCO (Pan et al., 20006) have revealed enhanced density of 
states at energies near ±7 mcV (i.e. well inside the bulk superconducting gap) 
with isotropic exponential dependence on distance from the core center. This might 
indicate that a field induced dx2_y2 +idxy order parameter is realized, see (Franz and 
Tesanovic, 1998; Pan ct al., 20006) and references therein. Tafuri and Kirtlcy (2000) 
have recently found spontaneous circulating currents around impurities in YBCO 
where the vortices contained fluxes of fractional size, in line with the idea of a complex 
order parameter induced near (magnetic) impurities (Graf ct al., 2000). Movshovich 
et al. (1998) found that the thermal conductivity of nickel-doped BSCCO suddenly 
dropped at about 200 mK, which was taken as evidence for a second superconducting 
transition to a state with dx2_y2 + id^y order. On the other hand, in a very recent 
experiment reported by Carmi et al. (2000), spontaneous magnetization of YBCO 
was found to appear at the superconducting transition temperature, meaning that the 
bulk order parameter may have complex (ij.2 _y2 +idxy symmetry. At this point, it seems 
crucial to collect more evidence for (or against) the dx2-y2 + idxy order and establish 
whether it appears in the bulk at Tc, in magnetic fields, or near imperfections like 
magnetic impurities at low temperature. 

Finally, the obvious, but perhaps naive, question concerns the all-important issue 
of the mechanism for high-Tc superconductivity, and the true symmetry of the order 
parameter. Despite the absence of a definitive theory of the pairing mechanism and 
symmetry of the order parameter - a consistent picture of d-wave singlet pairing and 
Andreev midgap states is slowly emerging, unifying a large body of experimental 
resiilts and providing a framework for critical examination of present and future 
experiments. This means, perhaps, that inconsistent results may be discarded and 
significant new effects discovered and accepted. If another paradigm shift becomes 
unavoidable, e.g. based on a consistent non-Fcrmi-liquid picture of the normal state, 
it is likely that this must incorporate the present notion of Andreev midgap states. 
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Appendix 

In this appendix we provide a simple derivation of the midgap states (MGS). In 
some detail, we first solve the Bogoliubov-de Gennes equation for bulk s-wave and 
rf-wave superconductors and then we discuss the clean NS interface and the concept 
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of Andreev reflection. The Andreev reflection phase shift used throughout the review 
is then obtained. Finally we show how the MGS are formed as a consequence of 
the sign difference between the lobes of the d-wave gap and why the formation of 
MGS is not sensitive to the detailed spatial dependence of the gap, although the MGS 
wavefunction itself depends on this shape. 

Some particularly useful references regarding the technique used to obtain the 
solutions of the Bogoliubov-de Gennes equation for different structures can be found 
in (Blonder et al., 1982; Bruder, 1990; Bagwell, 1992). For the s-wavc NS junctions 
we refer to Blonder et al. (1982) and for the d-wave case we refer to Bruder (1990). 
Bagwell (1992) derived the Andreev states present in s-wave SIS structures. 



Appendix A.l. Plane wave solutions of Bogoliubov-de Gennes equation 
The BdG equations for s-wave superconductors are 
r Hou{v) + A{v)v{v) = Eu{v) 

\ -Hov{r) + A*{r)u{r) = Ev{r), ^ ' 

where u and v are the electron and hole components of the wavefunction, Hq = 
— 7i^V^/2m — /X in the simplest case, and E is measured relative to the chemical 
potential /U. 

Equation (A.l) is easily solved for a bulk superconductor with A(r) = A = const 
by making the following ansatz: 

*s(r) = ( ) e^"--, (A.2) 

where u and v do not depend on r. The secular equation restricts the modulus of the 
k-vectors, k = nfc, to 



k = J^iii ± - A2) ^kF±^, (A.3) 



where ^ = s/E"^ — A^. In the last step we assumed that the quasiparticle energy 
is much smaller than the Fermi energy. There will then be two types of wave 
functions, electron-like and hole- like, depending of what type of solutions they describe 
in the normal limit (A — > 0). Labeling these by and k^ respectively, where 
k'^'^ = kp i a^/hvF (c = sgni?), we have 

^ tJO ^ ,,.= n.r^ ^ ^0^ ^..'"n-r^ ^^^^ 

where w,o and vq are given by vq/uq = (E — (7^)/A together with the normalization 
condition |woP + |vop = 1. Above we assumed that > A since there cannot be any 
propagating solutions for energies below the gap. However, there are also exponentially 
decaying or growing solutions which become important for surfaces. By defining 
^ = icrV A^ — E"^ for \E\ < A, the above expressions describe the exponentially 
decaying (in the direction of propagation) solutions for subgap energies. 

For the normal metal, the solutions are obtained from the above by letting A — > 0. 
The electron and hole wave functions can then be written as 

= ( J ) e''~"", = ( 1 ) (^•^) 

where k^'^ « fci? ± E/hvp- 
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Note that the quasiparticle's direction of propagation is given by the group 
velocity Vg = {l/h)dE/dk, which for a positive k is in the positive (negative) direction 
for electron-hke (hole-hke) quasiparticles. 

Appendix A. 2. d-wave case 

In the d-wave case, the Bogohubov-de Gennes equations are non-local 

Hou{r) + Jdr'Air,r')v{r') = Eu{r) 
~Hovir) + Jdr'A*{r,r')u{r') = Ev{v). ^ 

Making the quasiclassical approximation as described in e.g. reference (Bruder, 1990), 
we separate out the fast l/kp oscillations from the beginning by making the ansatz 

\ = ( ^(r) \ e^^— • (A.7) 

v{y) ) \ v{t) ) 

where the spatially slowly varying envelope functions u(r) and ■y(r) satisfy the 
quasiclassical BdG equation 

J — iTiWiTii • V{((r) + A(n, r)z;(r) = Eu{v) ( \ 9\ 

\ iTiWFn- V{i(r) + A*(n,r)u(r) = Ev{v) ^ > 

where terms of order l/kpS^^ and higher have been neglected. This approximation 
is valid when the inequality A <C Ep holds [see e.g. (Andreev, 1964)]. For the 
spatially independent order parameter, taking the r-dependence of u{r) and v{v) 
to be exp(iK(n)n • r), the secular equation restricts K(n) to ±^(11)/??,?;^?, where 
^(n) = ^E'^ — |A(n)p. In the same way as in the s-wave case, we then find electron- 
like and hole-like solutions 

n = ["i) e^'^-. *i = (;;;;) e^^^-. (a.9) 

where k'^'^ = kpi a^(n) /fivp ■ and w,o and wo are given by vq/uq = {E — cr^(n))/A(n) 
together with the normalization condition |moP + l^op = 1. 

As in the s-wave case, for siibgap energies we can find exponentially growing or 
decaying solutions given by the above expressions if we let ^ (n) = icr ■\/|A(n)P — E^. 

Appendix A. 3. Andreev reflection 

On a superconducting surface, electrons are reflected as holes and vice versa. In this 
section wc will discuss this basic process more quantitatively by solving the BdG 
equation for the normal metal/s-wave superconductor interface. Putting the N/S 
interface at a; = and assuming that the gap depends on the space coordinate as 

- { Ae-, 1 1 1 (^-^O) 

we are left with the simple problem of matching elementary plane wave solutions of 
the BdG equation in each region at the surface. The plane wave solutions of the BdG 
equation are derived above in this appendix. Assuming, for simplicity, propagation 
along the x-axis only, we write down an ansatz wavefunction describing an electron 
incident on the superconductor from the normal metal side 

= ( M e^fc«^ + ^ J J J e^^"^ + B ( J j e-''^"^, (A.ll) 
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where we have taken into account that the electron may be reflected as an electron as 
well as a hole. On the superconductor side we include transmitted electron-like and 
hole-like quasiparticles 

*»=^( j-<0 "c^-'O 

At the interface (at a; = 0) the wave function continuity gives 

B = D = 0, C = 1/uo, A = {vQ/uo)e-'^. (A.13) 

During this process of Andreev reflection, the electron incident from the normal metal 
is reflected as a hole with the probability amplitude A. The probability for Andreev 
reflection is |Ap = 1 for energies within the gap while quickly decreasing outside the 
gap. It is then convenient to write the amplitude of Andreev reflection for \E\ < A 
on the form 

_. _. E 
A = ae ^■^e 7 = arccos— , (A.14) 

which defines the Andreev reflection phase shift 7 used in the text. Note that also 
the phase x of A is picked up during Andreev reflection. Charge conservation during 
Andreev reflection is preserved by letting a Cooper pair enter the superconductor. 
Consequently, a charge 2e is transferred from the normal metal to the superconductor, 
which explains how subgap current may appear. In the same way, a hole incident 
on the superconductor from the normal metal side may be Andreev reflected as an 
electron transferring a charge 2e from the superconductor to the normal metal. The 
phase shift is in this case — 7 + x as can be checked by repeating the above calculation 
for an incident hole, keeping the rest of the ansatz as it is written. 



Appendix A. 4- Midgap states 

In this section we discuss quasiparticle reflection from the insulator/d-wave 
superconductor interface for energies within the superconducting gap and the 
formation of MGS. In the text we discussed MGS having the structure in figure 3 in 
mind. Here we show an alternative way of finding the MGS, in which it becomes clear 
that the normal metal plays no crucial role: the MGS is formed only as a consequence 
of the sign difference between the d-wavc lobes. When the size of the normal metal 
in figure 3 becomes vanishingly small, the remaining part of the MGS wavefunction is 
the exponential tails in the superconductor. Here we study this wavefunction first for 
step-function shape of A and then for arbitrary shape. 

Specular refiection at the surface imposes the following boundary condition for 
functions in equation (A. 9): the eigenstate consists of a superposition of states 
with trajectories n and n, and the total wave function is zero at the surface, i.e. 

*s(n,0) = -*s(n,0). (A.15) 

For spatially constant A, the solutions of equation (A. 8) with E = and decaying at 
infinity are given by 

*s(n,r) = ( ^ \e^kFn■r-\A/hv,r^\x^ 

\isgnA / 
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If A = — A, these functions have similar vector structure and they can be matched at 
the surface, a; = 0, according to equation (A. 15), and therefore there are MGS in this 
case. 

Equation (A. 8) has solutions at E = also for arbitrary spatial variation of the 
gaps A(r) and A(r). By introducing the following new variable. 



c 



: j dl\A{l)\/nvF, (A.17) 

where I = n r, one can reduce equation (A. 8) to: 

du 

-i— + sgn{A)v = 

dv^ (^-18) 
«— +sgn(A)u = 0. 

These equations have solutions similar to the ones in equation (A. 16), which leads to 
MGS. 
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Figure 1. Quasiclassical path giving the bound states at a specular surface of 
a superconductor with anisotropic s-wave order parameter A (6) (schematically 
drawn in the right part of the figure). N is a normal region of size L ~ to 
model the gap suppression near the surface. The quasipaxticles are trapped in the 
surface region because of normal reflection at the vacuum/metal (IN) surface and 
subgap Andreev reflection at the NS surface. The angle 9 gives the propagation 
direction relative to the surface normal. During Andreev reflection, an electron 
(e) is converted into a hole (h) with the phase gain —7 — x- According to the 
Bohr-Sommerfeld quantization condition, the total phase accumulated during one 
cycle is equal to 2nn. 
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Figure 2. (a) Density of surface states at an s-wave superconductor surface with 
a normal metal overlayer of thickness 1.5^0. The peak energies are the bound 
states for trajectory angle = and the tails towards lower energies are due to 
the angle integration; for increasing angles the energy of the surface states are 
continuously shifted down relative to the 9 = state, (b) The same for a gap of 
d-wave synmictry. The crystal is rotated by 1.5° relative to the surface normal so 
that midgap states are formed (the delta peak at zero energy). 
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Figure 3. Quasiclassical path giving the bound states at a specular surface 
of a two-dimensional d-wave superconductor with order parameter A(9) = 
Aocos[2(6 — a)] where is the propagation angle, compare Fig. 1, and a is the 
orientation angle of the order parameter relative to the surfax;e normal. N is a 
normal region of size L ~ to model the gap suppression near the surface. 
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Figure 4. Temperature dependence of the penetration depth A(oc S~^) measured 
on YBCO films containing differently oriented linear defects introduced via ion 
bombardment. Due to the presence of MGS at the defects, which respond 
paramagnetically to an external field, there is an upturn of A at low temperature. 
The upturn in the unirradiated film is thought to be due to MGS at internal 
interfaces (grain boundaries). What is important is the enhancement of the 
paramagnetic effect with increasing orientation angle, compare the angle a in 
Fig. 3, from to 45° of the defects. Prom Walter et al. (1998), copyright (1998) 
by the American Physical Society. 
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Figure 5. Quasiclassical trajectory illustrating the Andreev bound state in 
a specular two-dimensional superconductor junction. Sl and Sft arc the left 
and right ci-wave superconductor electrodes with orientation angles aj, and an 
respectively. The junction region consists of a barrier (insulator I) and normal 
regions Wj, and Nj^ which may be due to the gap suppressions. The Andreev 
state can be thought of as a hybridization of the surface states (Fig. 3) of the left 
and right sides of the insulator. 
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Figure 6. (a) Dispersion relation E{<p) = A(6) D{9) cos for the Andreev 
bound states in the speculaj: d-wave junction with orientation d^/i/d^/i for the 
injection angle 6 = 7r/4. The transparency is D{6 = 7r/4) = 0.01. (b) Angle 
integrated current-phase relation for the same junction for very low temperature 

(T = 0.003Tc). (c)-(d) The same but for the orientation d^ In this case 

the dispersion relation is E{ij>) = A(S)-^/D(^sin()i and the current is 7r-shifted 
relative to the ci7r/4/<i7r/4 case. The angle averaged transparency is < 13 >= 0.026. 
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Figure 7. Critical current as a function of temperature for tlircc different 
orientations of the d-wave superconductors: d^/4/i^,r/4 (solid line), a!,r/4/'^-7r/4 
(dashed line), and d^/g/d_^/g (dot-dashed line). The negative sign of Ic indicates 
that the junction minimum is at </> = tt (tt junction). The — > tt junction 
crossover with decreasing temperature for the d^^g/d_^/g junction is due to 
the competition between the MGS currents (dominating at low temperature) 
and current contributions from gap-edge bound states and continuum states 
(dominating at high temperature). The angle averaged junction transparency 
is <D >= 0.026. 
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Figure 8. Critical current for the d,r/4/'^7r/4 orientation. At low temperature, 
Ts, a subdominant component of the order parameter with s-wave symmetry is 
assumed to appear and be 7r/2 out of phase with respect to the dominant (i-wave 
component. This results in a cut off of the MGS contribution and the junction 
changes from a 7r-junction to 0-junction when the temperature is decreased. 
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Figure 9. (a) Andrccv bound states in a specular d-wave junction witli tlie 
a!o/(i^/4 orientation for tlie injection angle d = it/9 and transparency Did = 
7r/9) ~ 0.01. The solid and dashed lines are the +ky and —ky states respectively, 
(b) The angle integrated current-phase relation for this orientation at three 
different temperatures: T = {0.001, 0.01, 0.1}Tc. (c) Phase dependence of the 
surface current density calculated to the right of the barrier for the same three 
temperatures. As the temperature is decreased, the ifcy MGS become unequally 
populated and start to dominate the Josephson current which then is enhanced 
and is of order D. At the same time the surface current becomes appreciable. 




Figure 10. Ic{T) dependences for high transparency, (D) = 0.68 [Z = 0.5), 
for three orientations: dg/do dashed line, <i,r/4/'^— 7r/4 dash-dotted line, and 
d—iso jdzQo solid line. The last orientation corresponds to the nominal orientation 
in the experiment by Arie et al. (2000). 
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Figure 11. Illustration of (a) the single-particle (quasiparticle) and (b) the two- 
particle (Andreev) processes in the NIS junction. At zero temperature allowed 

processes connect occupied states below zero energy (the Fermi level) in the 
normal metal with (a) unoccupied states above the gap in the superconductor, and 
(b) above zero energy in the normal metal. We have aligned the Fermi energies of 
the two sides by a gauge transformation in the right superconductor; consequently, 
quasipajrticles are accelerated by eV at each passing through the junction region. 
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Figure 12. Current voltage relations and corresponding differential conductances 
for the N/da junction with two different orientations: a = 0, shown in (a)-(b), 
and a = 7r/4, shown in (c)-(d). The zero-bias anomaly in (c)-(d) is due to the 
MGS resonance. The barrier is modeled by a (5-function potential with a large 
tunnel cone with angle averaged transparency < D >= 0.026 which emphasize 
the importance of the MGS resonance which appear for finite injection angles $ 
only. The temperature is zero and ctq is the normal state conductance. 
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Figure 13. Conductance spectra, normalized by the normal background 
conductance, obtained on YBCO at 4.2 K by scanning tunneling microscopy 
fitted by the scattering theory of MGS, assuming the a = 7r/4 orientation. The 
parameters of the theory are the maximum d-wave gap Aq = 17.8 meV, the 
normalized barrier height Z = 2mH/h^kp = 3 (H is the height of the barrier 
potential), and tunnel cone parameter = 10. From Kashiwaya et al. (1995), 
copyright (1995) by the American Physical Society. 
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Figure 14. Demonstration of the doppler shift of MGS in the presence of a 
magnetic field H = {0.2, 0.4, 0.6, 0.8, 1.0, 2.2, 3.4, 5.0}T in YBCO/Pb junctions, 
where the curve at the top is for the lowest field. The spectra are normalized by 
the high-voltage conductance at H = 0.2T and fitted by a Lorentzian centered at 
the bound states energy Eg = vppg cosff averaged over angles 6 with a tunnel 
cone D(8) = 1 for \9\ < 8c = 15° and D(9) = otherwise. The fit parameters are 
the Lorentzian width F = [Ao/A(9)] X 0.8 meV and vpps- In inset (a) theoretical 
results are shown for a larger tunneling cone for comparison. In inset (b) the 
maximum Doppler shift Em = vpPs cos(7r/2 — 6c) and the ZBCP splitting 5 are 
both shown as a function of applied field H. Prom Aprili et al. (1999), copyright 
(1999) by the American Physical Society. 
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Figure 15. (a) Conductance of YBCO/Ag ramp-edge junctions showing the 
orientation dependence of the zero-bias conductance peak as predicted by the 
MGS scenario: for a = the peak is absent [the (100) curve], while it is present 
for the a = 7r/4 orientation [the (110) curve], (b) Temperature dependence of 
the conductance peak [the (110) curve in (a)]. The peaJc height scales as 1/T, 
as shown in the inset, while the width is temperature independent. Prom Wang 
et al. (1999), copyright (1999) by the American Physical Society. 
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Figure 16. Conductance spectra of Zn-doped BSCCO obtained by scanning 
tunneling microscopy in the c-axis direction at 4.2 K. The spectrum taJsen at a Zn 
impurity site contains a huge zero-bias anomaly (open circles) , while the spectrum 
taken at a clean region is of the usual d-wave form with peaJfs at the gap voltages 
(arrows). Reprinted by permission from Nature [Pan et al. (2000a)] copyright 
(2000) Macmillan Magazines Ltd. 
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Figure 18. Illustration of the (a) single-particle, (b) two-particle, and (c) three- 
paxticle processes in the SIS junction. At zero temperature allowed processes 
connect occupied states below the gap in the left superconductor with unoccupied 
states above the gap in the superconductors. We have aligned the Fermi energies of 
the two sides by a gauge transformation in the right superconductor; consequently, 
quasipajrticles are accelerated by eV at each passing through the junction region. 
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Figure 19. The total dc current in junetions with orientations (a) da/di), (b) 
^/d-n/A where As = 0.2A() (As is the s-save superconducing gap and Ao is the 
maximum d-wave gap), (c) do/d^n, and (d) d^n/d^n. The transparency of the 
tunnel barrier is < D >= 0.026 in all cases and we assume zero temperature. 
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Figure 20. (a)-(c) The total current [the same curve as in figure 19(b)], the 
one-particle current, and the two-particle current in the sfd^jn junction with 
Ag = 0.2Ao where A^, is the s-save superconducing gap and Aq is the maximum 
d-wave gap. The peak at the gap voltage is due to the MGS resonance on the right 
side overlapping with the gap-edge resonances at the entrance and exit points on 
the s-wavc superconductors side, sec inset in (c) where the circles indicate presence 
of resonances. Equivalently, we can say that the large BCS gap-edge density of 
states pushes more current through the MGS, compared to the NS ease, (d)- 
(f) The total current, the one-particle current, and the two-particle current in 
the I ^1 I ^ junction. The onset at the gap voltage is clue to the hare MGS 
resonance, see inset in (f) where the circle indicate presence of resonance. The 
transparency of the tunnel barrier is < D >= 0.026 in all cases and we assume 
zero temperature. 



